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Abstract: We propose a picture for the UV properties of Galileon field theories. We con¬ 
jecture that Galileons, and all theories incorporating the Vainshtein mechanism, fall into 
Jaffe’s class of ‘non-localizable’ held theories characterized by an exponential growth in 
their Kallen-Lehmann spectral densities. Similar properties have been argued to arise for 
Little String Theories and M-theory. For such theories, the notion of micro-causality and 
the time ordering used to dehne the S-matrix and correlation functions must be modihed, 
and we give a Lorentz invariant prescription for how this can be achieved. In common with 
General Relativity (GR), the scattering amplitudes for Galileons are no longer expected to 
satisfy polynomial boundedness away from the forward scattering or hxed physical momen¬ 
tum transfer limits. This is a rehection of the fact that these theories are fundamentally 
gravitational and not local held theories. We attribute this to the existence of a local¬ 
ity bound for Galileons, analogous to the Giddings-Lippert locality bound for GR. We 
utilize the recently developed Galileon duality to dehne a UV hnite, Lorentz invariant, 
quantization of a specihc Galileon theory for which the energy of all states are positive def¬ 
inite. We perform an explicit computation of the Wightman functions for this theory, and 
demonstrate the exponential growth associated with the locality bound. In analogy with 
GR, the bound is correlated with the absence of Galileon Duality (i.e. Diffeomorphism) 
invariant local observables. We argue that these theories can be quantized in a manner 
which preserves Lorentz invariance and macro-causality and that the latter ensures that 
the superluminalities found in the low energy effective theory are absent in the full theory. 
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1 Introduction 


The discovery of cosmic acceleration, and the desire to develop theories of gravity in 
which the dynamics is distinct from general relativity in the infrared in order to tackle 
this problem, has led to the development of new classes of field theories incorporating 
screening mechanisms. A whole class of snch theories incorporate a gravitational screen¬ 
ing mechanism, known as the Vainshtein mechanism [1], whose existence is necessary for 
the observational viability of these theories. However, the status of these theories as con¬ 
sistent quantum field theories has remained uncertain. No explicit UV complete model 
which incorporates the Vainshtein mechanism is known, and there have been a number 
of arguments to suggest that such a completion may be impossible. The simplest model 
incorporating the Vainshtein mechanism is what is now known as the cubic Galileon model 
[2]. It arose as the effective description of the helicity-zero mode of the graviton in the 
Dvali-Gabadadze-Porrati (DGP) model [3]. Specifically it was defined as the graviton mass 
m —)• 0 double scaling/decoupling limit of the DGP model in which the dynamics of the 
helicity-two modes decouple, Mpjanck oo, whereas the low energy gravitational dynam¬ 
ics of the helicity-zero mode vr survive. In this limit an accidental non-linearly realized 
Galileon symmetry tt —)■ tt + c -h survives as a remnant of the higher dimensional 
diffeomorphism symmetry. Related to this, vr may be viewed as a Goldstone model for a 
spontaneously broken higher dimensional spacetime symmetry group (see Sec. 5.3). 

The authors of [4] then generalized this structure to the most general scalar theory 
consistent with the same non-linearly realized Galileon symmetry which was intended as an 
effective description of some unknown infrared completion of gravity. It was subsequently 
realised that the entire structure of Galileon theories arise naturally as decoupling/scaling 
limits of Massive Gravity theories [5], i.e. ‘hard mass’ Massive Gravity theories (like ‘soft 
mass’ DGP) are examples of IR completions to Galileons, and this structure was impor¬ 
tant for determining the unique ghost-free theories of massive gravity [6]. The Galileon 
structure has since been shown to arise in a whole host of infrared modifications of gravity 
which rely on giving a hard or soft mass to the graviton in a Lorentz invariant way (see 
[7] for a review on some of these approaches). More recently, the Galileon has been shown 
to play a special role in soft limits of scattering amplitudes [8], and the full tree-level scat¬ 
tering amplitudes of a special Galileon have been given in [9] using a novel ‘dimensional 
reduction’ from d + d dimensional Einstein gravity. This special Galileon has been shown 
to admit an additional symmetry (beyond the Galileon symmetry) in [10]. 

All of these developments have taken place at the level of the Low Energy Effective 
Eield Theory (LEEET) (or tree-level S-matrix [9]) with the understanding that some UV 
completion exists. However at the same time, by treating the Galileon decoupling limit as 
an ordinary field theory, doubt has been cast on the existence of any such UV completion 
[11]. In a standard field theory, the key properties are that it exhibits a UV completion 
whose S-matrix is Unitary, Analytic, Lorentz Invariant and Polynomially Bounded (in the 
form of the Eroissart bound). Indeed any field theory satisfying the Wightman axioms will 
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automatically satisfy polynomial boundedness and analyticity of the S-matrix, at least for 
fixed physical momentum transfer [15, 64, 65]. It is easily shown that the cubic Galileon vi¬ 
olates one (or more) of these properties [11] and that same argument extends to all Galileon 
models, and by extension Massive Gravity [6] and its Multi-gravity extensions [12, 13]. 

We will argue here and elsewhere [14] that it is the polynomial boundedness assumption 
that is not appropriate for Galileons^. Our assumption will then be that the Wightman 
functions are no longer tempered distributions, which amounts to the UV description of 
Galileons admitting some degree of non-locality. In fact it is generally believed that Gen¬ 
eral Relativity and string theory violate polynomial boundedness, at least for fixed angle 
scattering or unphysical momentum transfer, in the region of small impact parameters for 
which black hole production can take place. This gravitational non-locality is intimately 
connected with the production of black holes controlling the high energy properties of 
scattering amplitudes and with the absence of off-shell local observables in a gravitational 
theory. We will argue that Galileon theories are in this sense fundamentally gravitational, 
with the Galileon duality playing the same role as diffeomorphism invariance in forbidding 
local off-shell observables. 

The idea that a theory whose infrared limit does not contain any dynamical metric 
and is a local effective field theory, should nevertheless be viewed as a gravitational theory, 
in the sense that it exhibits some high energy non-locality, is unusual, but is not new. An¬ 
other well known example are Little String Theories (LSTs). Despite being field theories 
in the infrared limit, their UV description will not be a field theory with a UV fixed point. 
We will elucidate these points below. In [16] it is argued that if their exists a field theory 
description of LSTs, then they are Jaffe type field theories of quasi-local type [17, 18]. Here 
we will argue that analogous properties are true for Galileons, i.e. that the UV Galileon 
fields are Jaffe fields of non-localizable type [19], with the only distinction from LSTs being 
the degree of non-localizability as defined by the order of growth of the Kallen-Lehmann 
spectral density. A related example, which we argue has many properties similar to out 
conjectured UV description of Galileons, is the ‘asymptotic fragility’ proposal for the de¬ 
scription of infinitely long strings [20, 21]. 

In short, we conjecture that the UV description of Galileons, and by extension all the¬ 
ories of Massive Gravity, Bi-gravity and Multi-gravity violates the condition of polynomial 
boundedness, but they respect Unitary, Lorentz Invariance and crucially Analyticity (the 
limitations imposed by the latter will be discussed elsewhere [14]). The main arguments of 
this paper are that: 

• We give evidence for the violation of polynomial boundedness by computing the exact 

^Polynomial boundedness is the assumption that the growth of the scattering amplitude as a function 
complex momenta ki is bounded by a polynomial, i.e. A{{ki]) < C| Polynomial boundedness 

follows automatically for theories in which the Wightman functions are tempered distributions, however it 
also follows for all strictly localizable fields [15] (at least for fixed physical momentum transfer) for which 
the Wightman functions grow slower than a linear exponential W{{ki}) < 
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Kallen-Lehmann spectral density of a specific Galileon and demonstrating its expo¬ 
nential growth. To do this we use the Galileon duality to define a specific quantization 
of a Galileon theory which is dual to a free theory. This fully quantum equivalence 
between non-localizable field theory and a localizable theory, which extends to the 
full S-matrix [14] is an example of the Borchers equivalence, for which the S-matrix 

and asymptotic states are identical [22]. For a Galileon in d dimensions we find 
p{E) ~ g(E/A)('*+2)/('i+i)_ 

• The exponential growth of the spectral density is consistent with semi-classical argu¬ 
ments for which the classical action scales as |5| demonstrating 

that much of the high energy behavior can be anticipated from semi-classical consid¬ 
erations [23-26]. It arises because the spectral density is dominated by intermediate 
states with a finite number of particles N for which N rsj (<^+2)/(li-i-i) _ 

• As a consequence of the exponential growth in the spectral density, we argue that 
Galileons fall into Jaffe’s class of non-localizable field theories. We further conjecture 
that all field theories exhibiting the Vainshtein mechanism are non-localizable held 
theories. 

• We argue that Galileons exhibit a locality bound in space and time which is a property 
of all non-localizable held theories. This is connected with the fact that correlation 
functions are no longer tempered distributions and must be dehned using a space of 
test functions which includes no functions of compact support. 

• Although Galileons break strict micro-causality/locality, we argue that they never¬ 
theless preserve macro-causality/locality. We give a prescription that dehnes macro- 
causality for the S-matrix and two-point correlation function out of vacuum, and 
argue that this implies the absence of superluminal propagation seen in the LEEFT 
in the UV theory. 

• It is known that non-localizable theories can exhibit standard properties, such as 
cluster decomposition, GPT, and an LSZ formalism [19], but are non-Wilsonian in the 
sense that they exhibit no UV fixed point, and no local operator product expansion. 

• We give a Lorentz invariant prescription for defining time ordered correlation func¬ 
tions and generalized retarded products in terms of Wightman functions. This pre¬ 
scription makes use of entire analytic functions of momenta having the same order 
of growth at large complex momenta as the Wightman functions. 

• The absence of off-shell, local observables is consistent with the viewpoint that 
Galileon theories are fundamentally gravitational. This is consistent with their ori¬ 
gin as decoupling limits of resonance (soft) or hard mass gravity theories, and the 
existence of the Galileon duality transformation which is simply a remnant of diffeo- 
morphism symmetry that survives in the decoupling limit. 
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• Our picture of the UV properties of Galileons is consistent with the ‘asymptotic 
fragility’ proposal [20, 21], in the sense there is no UV fixed point, and the ‘clas- 
sicalization’ proposal, in that sense that the high energy scattering is expected to 
be semi-classical [23-26] and scattering should be dominated by the production of 
V-particle states with N ~ 


This proposed picture for the UV properties of Galileons is radically different than 
standard low energy effective field theory expectations based on the idea that the Galileon 
strong coupling scale A is the cutoff of the effective field theory (EFT). In particular, no 
new particle states need to arise at the scale A in order to resolve unitarity. Rather, A is 
a fundamental scale, intrinsic to the definition of the UV completion and the properties 
of the non-localizable fields^, a role played analogously to the Planck scale in quantum 
gravity. There is no UV renormalization group fixed point, i.e. no conformal behavior at 
high energies. 

Nevertheless, at energies E <C A, the conjectured UV completion will be consistently 
described by the Galileon LEEFT. In a perturbative expansion in E/A, the exponen¬ 
tial growth becomes polynomial, thus recovering local physics at low energies and large 
distances. Specifically as long as we look at separations jxj S> r^{E), where r^{E) = 
is the Vainshtein/classicalization radius, and E < A the LEEFT will be 
satisfied. As soon as E > A, we have r,^[E) > 1/A and there will be a region even at 
distances > 1/A where the LEEFT locality properties will be violated and there will be 
some manifestation of the non-localizability. These properties are analogous to the those 
expected for GR, as we discuss in Sec. 3. 

2 Galileons as Effective Field Theories and Beyond 

The traditional effective field theory (EFT) point of view on Galileons can be summarized 
by the usual EFT paradigm: identify the low energy degrees of freedom and write down 
every local operator consistent with the low energy symmetries. In the case of Galileons the 
symmetries are linearly realized Poincare invariance and the non-linearly realized Galileon 
and shift symmetries tt —?• vr -|- -|- c. These symmetries may be formalized via the 

coset construction (see Sec. 5.3) but the simplicity of the Galileon symmetry means that 
it is easy to write down the form of the most general interactions without this construction. 


^We shall continue to nse the historical terminology ‘non-localizable’ rather than ‘non-local’ to dis¬ 
tinguish these models from other distinct non-local field theories fonnd in the literature. For instance a 
common approach in non-local field theories is to modify the Feynman propagator by an entire function 
without introducing an additional pole with the hope of improving Euclidean perturbation theory. From 
our perspective such a modification violates unitarity since the associated Wightman functions do not sat¬ 
isfy the positivity requirement of the Kallen-Lehmann spectral representation, whereas all the Jaffe type 
non-localizable theories respect unitarity in the form of a positive Kallen-Lehmann spectral representation. 
Furthermore we make no claims for improved perturbation theory, here the non-localizability only arises in 
the non-perturbative regime. 
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There are two classes of interactions that arise for the Galileon, those built out of the 
manifestly Galileon invariant combination j and a finite number of special 

Galileon interactions whose Lagrangian transforms as a total derivative under the Galileon 
transformation. Thus we may write the Wilsonian effective action for the Galileon as 

5'Galileon = S'a + (2.1) 

where Sa are the lower derivative ‘Wess-Zumino’ combinations (index contraction sup¬ 
pressed) 

/ di 

d-^xAV^a^eenV" (2-2) 

n=0 

{a = d/2 + 1), € is the Levi-Civita tensor, and 5 b contains the infinite number of manifestly 
Galileon invariant combinations which we write schematically as 

/ a2n 

d'^xA'^+i , (2.3) 

n,m 

and a term of the form includes all scalar contractions of all combinations of 

derivatives applied to different II’s. There exists a non-renormalization theorem that vr 
loops to not renormalize Sa [27, 28], however the terms Sb are renormalized. The non¬ 
renormalization theorem implies that the strong coupling scale A and coefficients do 
not flow, but it does not improve the regime of validity of the perturbative LEEFT. In a 
Wilsonian cutoff approach, all terms (3n{-^w/^) rnust be included and will depend on the 
Wilsonian cutoff Aw 

As is well known, the special finite number of contributions 5 a lead to second order 
equations of motion. For this reason cosmological applications of Galileon theories usually 
focus on these finite number of terms. However, there is no requirement in an EFT that 
we should truncate to only those operators that have second order equations of motion. 
Rather all operators consistent with the symmetries should be included, and the would be 
Ostrogradski ghosts associated with the higher derivatives will have masses at what would 
then be the cutoff of the EFT Acutoff ~ A. Thus they lead to no unitarity violation in the 
regime of validity of the EFT and all higher derivative operators can be dealt with in a 
perturbative sense. 

Assuming all the coefficients an and (3n,m are of order unity, then the regime of validity 
of the EFT theory may be taken to be 


5 <C A (2.4) 

n<l. (2.5) 

We do not require vr <C A because of the Galileon and shift symmetries. To be precise as 
long as the shift symmetry is exact or only approximately broken, a large vev for tt will 
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not generate large quantum corrections as long as the above conditions (2.4) are met. 


This regime may be qualitatively understood as the regime in which we expect the 
expansion of the effective action to converge. In reality, all EFT expansions are asymptotic 
and so have technically zero range of convergence, however this regime is the one in which 
we may reliably truncate the expansion to the first few terms. 

2.1 Pushing the boundary of the EFT 

Although the previous discussion accounts for the regime allowed by a traditional effective 
field theorist, it has been argued that we may be able to push beyond this scale. A simple 
counting of powers shows that all the terms in S'b have at least one more power of derivatives 
than Sa (see for example [29]), i.e. 

Sb ~ ^5'a • (2.6) 

Thus as long as 5 <C A it is conceivable that we may focus only on the Sa part of the 
action, at least for determining semi-classical solutions. It may be plausible that it is even 
possible to have classically ft ~ 1 with the derivatives providing the expansion parameter 
necessary to make sense of the EFT. This is a rather modest improvement of the range, 
but does not lead to any significant changes in the phenomenology. 

A more radical proposal goes under the umbrella of the Vainshtein mechanism [1], 
and is central to the phenomenological viability of Galileon and massive gravity theories. 
Here the idea is that we may plausibly take backgrounds for which H ^ 1 and still remain 
in the regime of validity of the EFT^. The argument for this is based on the following 
reasoning: First take the effective action to include only those terms from Sa- Then 
determine a classical background solution which follows from the equations of motion for Sa 
which allows for a Vainshtein or strong coupling region for which H ^ 1. This background 
generically renormalizes the kinetic term for fluctuations in such a way that the effective 
strong coupling scale becomes redressed to something which is parametrically larger than 
A 

A* ~ Z^A , (2.7) 

where /3 is a positive model dependent parameter and Z ^ 1 is the effective wave function 
renormalization [30-32]. In this way we may argue that around a background with H S> 1 
that the cutoff of the effective field theory is really A* and so we may now push into the 
new regime of validity 

5 <C A* . (2.8) 

If Sir denotes the fluctuations around the background solution, then we may similarly argue 
that we can trust the reorganized EFT in the regime 

dddir <C A* . (2.9) 

®For this it is essential to distinguish the strong coupling scale A from the cutoff of the EFT. 
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2.2 Problems with the EFT description 

The Vainshtein proposal requires a very special behaviour of the UV completion of the 
Galileon theory. In a generic EFT, we would expect the scale A to be associated with the 
mass of some new degree of freedom. If this new degree of freedom exists at this scale 
then we would expect that unitarity would require us to include it if we begin looking at 
energies up to the scale A* > A. The Vainshtein mechanism implicitly assumes that this 
new degree of freedom does not need to be included which is generically false. 

To see this more concretely we just need to do the previous calculation in a different 
order. Suppose we begin with the total action for the Galileon including S'b. We then look 
for a background solution of this action which exhibits 11 ^ 1 and perturb around it. Even 
after we have perturbed around this solution there will exist the following infinite set of 
contributions to the effective action which are schematically 

d'^xA'^+i . (2.10) 

n 

Even with the wave function renormalization, the expansion of these infinite set of terms 
will breakdown at the scale 9 ~ A. In fact it is precisely these terms which would be 
associated with a modification to the Galileon propagator due to the existence of other 
massive poles. The point is that wave function renormalization may change the scale of 
interactions but they do not change the position of poles in the propagator. A resolution 
to this is to imagine that the UV completion is such that /3n,2 all happen to vanish (or be 
signihcantly smaller in magnitude). This is conceivable, but implies an extremely nontriv¬ 
ial requirement on the UV theory for it to be consistent with the Vainshtein mechanism at 
low energies. 

A perhaps more serious problem to the validity of the Galileon EFT description both 
at strong and weak coupling is the fact that it can be shown that no theory whose UV 
completion satisfies Lorentz invariance, unitarity, locality (polynomial boundedness) and 
causality (analyticity) can give rise to a Galileon EFT at low energies [11].^ Since every 
local quantum field theory satisfies these requirements there can be no local field theory for 
which Galileons are the low energy EFT, which means they are arguably not of interest as 
EFTs! It is conceivable of course that they admit a Lorentz violating UV completion (for 
suggestions along these lines for massive gravity see [33]), but then there was no reason 
to write down a Lorentz invariant theory in the first place. In fact the space of Lorentz 
violating Galileon theories is much larger due to the weakening of the symmetries. As a 
result of these arguments many people have viewed Galileon theories, and by extension 
massive gravity theories, as sick from the get go. It is our intention here to show that such 
a conclusion would be too quick, and results from interpreting Galileon type theories as 
local field theories, rather than intrinsically gravitational theories. 

^The original arguments were made for the cubic Galileon model, but it may be easily show that these 
apply to the generic Galileon models. 


- 8 - 



2.3 Beyond Effective Field Theories? 

Although not immediately apparent, the two previous problems are connected. It may be 
shown that in every theory known so far which exhibits classically the Vainshtein mech¬ 
anism, the S-matrix is such that it does not admit a standard ‘Wilsonian’ or ‘local field 
theory’ UV completion at the scale A. Thus the Vainshtein mechanism appears to be 
incompatible with a standard assumption that the scale A is associated with new states 
in a local field theory UV completion. This tells us that in these classes of theories we 
should view the UV completion in different terms [34]. It will be our contention here that 
Galileons can be understood as a class of field theories considered in the 1960/1970’s which 
are known as ‘non-localizable’ field theories. These field theories are not Wilsonian in the 
sense that they are not dehned via UV RG fixed points, and their operators do not sat¬ 
isfy a standard operator product expansion, they may nevertheless satisfy other standard 
properties such as cluster decomposition, analyticity, LSZ construction, crossing symmetry 
and, CRT symmetry [19]. Before getting to a discussion of this perspective we first clarify a 
number of arguments that have been made pertaining to the UV properties of these theories. 

One resolution to the problems of the EFT of Galileons is to imagine that view¬ 
ing Galileons as EFT theories, whilst not technically inconsistent in the EFT safe region 
9 <C A, n <C 1, is nevertheless fundamentally misleading since it assumes a standard UV 
completion at the scale A which is known not to exist. Thus a more radical proposal, but 
one which is consistent with the Vainshtein mechanism, is that A is not the cutoff of the 
Galileon theory, but a strong coupling scale above which the theory self-unit arizes without 
the introduction of new degrees of freedom^. This is the essence of the ‘UV completion 
via Glassicalization’ idea of [23-26] and we shall see that many aspects of this picture are 
consistent with the properties of non-localizable field theories. To understand these ideas 
it is worth commenting on what we mean by the breakdown of the EFT. 

Let us for now focus on the hnite Galileon terms S'a and ask the question, at what 
scale does unitarity breakdown? There is a well-known criterion for this, tree-level unitarity 
breaking, which proceeds as follows. It is usual to decompose the S-matrix in the form S = 
1 + iT and then compute the momentum conservation stripped transfer matrix {f\T\i) = 
{2'K)'^5^{pf — Pi)A{i —>■ /). Unitary is encoded in the optical theorem 

i(r+-r) = rtr ( 2 . 11 ) 

which between states amounts to the generalized optical theorem schematically given as 
f)* - Mf *)) = - Pi)A{i n)A{n f) . (2.12) 

n 

These equations impose bounds on the magnitude of the amplitude. Specifically, for 2 —)■ 2 
scattering, expressing the amplitude in partial waves 

A{s,t{cose)) = J + l)-Pi(cos6l) —(az(s) - 1), (2.13) 

V z=o * 

®By self-unitarization we mean that there are no new poles in the propagator. 
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then unitarity amounts to the condition that |az(s)| < 1. 

When computing the scattering amplitude A(s, t, u) for 2 —)• 2 elastic scattering to 
tree-level for a Galileon, we would find in terms of the standard Mandelstam variables, 
schematically ^(s, t, u) ~ (s^ -|- Expressing this in terms of partial waves this 

may easily be shown to violate the optical theorem when ^/s ~ A (assuming A ;§> m)®. 
For this reason it is usual to regard A as the cutoff of the Galileon EFT, meaning that we 
cannot reliably trust the EFT at and above this scale since without new physics coming to 
the rescue it seems sure to violate unitary. 

However this violation of the optical theorem is misleading for the following reason. 
We can equally well in perturbation theory express the S-matrix in terms of the reac- 
tion/reactance/Heitler matrix K as 


^ _ 1 - iK/2 

l + iK/2 

Similarly K may be expressed in a momentum conservation stripped form 

'\f\K\i) = - Pi)K.{i /). ( 2 . 15 ) 

In this language, unitarity and the optical theorem is simply the statement that K is 
Hermitian which amounts to the linear condition 

JC{i ^ f) = JCif ^ i)*. ( 2 . 16 ) 

If instead of computing perturbative corrections to T we compute them directly to AT, then 
at each order in perturbation theory we will find that K remains Hermitian b This follows 
automatically for any theory whose Hamiltonian is Hermitian. From this Hermitian K we 

®We note that although the Galilon symmetry forbids a mass term m, since this is only a global symmetry 
we expect a small mass to be present. In particular in the context of massive gravity the Galileon arises as 
a massive particle. Including a small mass is also important for making use of analyticity arguments since 
it guarantees a small window in the Mandelstam plane where the amplitude is analytic around the real axis 
which is needed to prove Hermitian analyticity [14]. 

^We may easily reorganize the standard perturbation theory expansion to directly compute corrections 
to K. For instance, in the interaction picture, with Hi the interaction, explicitly to second order in 
perturbations K is given by 

/ oo • n oo n oo 

dtHi(t)-- At At'e{t-t')[Hi{t),Hi{t')\ + ... . (2.17) 

-OO " J —oo J —oo 

This is reminiscent of the in-in formalism in that non-time ordered correlators arise through the com¬ 
mutator. This expression is manifestly Hermitian due to the properties of the commutator. The exact 
non-perturbative expression for K follows from the two equations [35] 

dtHi{t)Vit) , Vit) = 1 - X / dt'e{t - t')Hi{t')V{t') , (2.18) 

-OO " J —oo 

where e{t — t') = 9{t — t') — 9{t' — t). Indeed the largely forgotten goal of the Lippmann-Schwinger paper 
[35] was to give a variational technique to compute K non-perturbatively. 


(2.14) 
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may then reconstruct the unitary S. Thus unitarity is never violated at any finite order in 
perturbation theory. From the perspective of the Feynman diagram expansion this may be 
viewed as a resummation of an infinite number of diagrams corresponding to the expansion 

OO 

5 = (2.19) 

n=0 

Truncating this expansion would imply a loss of unitary at the same scale, but there is no 
need to truncate, nor is there any need to phrase this in terms of the Feynman diagram 
expansion. This resummation may be performed explicitly at the level of the partial wave 
amplitudes and is straightforward to do in the elastic limit, or almost elastic limit in which 
only a finite number of intermediate particles contribute®. 

An explicit example of this unitarization via resummation is the self-healing or self- 
unitarization mechanism of [37]. There it is shown that in QCD with arbitrary numbers of 
flavours Nf and colours Nc, it is possible to create a hierarchy using Nc and Nf such that 
tree-level unitarity breaks down for the low energy Chiral Perturbation Theory (ChPT) 
EFT at scales below the actual cutoff of the ChPT EFT, i.e. the scale of new physics 
a.k.a. QCD. This may be resolved by computing K (in its partial wave representation), 
combining to give S and showing that unitarity is maintained. At a practical level this 
corresponds to resumming the partial wave amplitudes into a manifestly unitary form. 
However there is physics in this resummation. In the K representation of the S-matrix it 
is clear that additional poles will generically arise already in perturbation theory due to 
the denominator. In the case of ChPT, the additional pole that arises is on the second 
Riemann sheet and so may be associated with a physical resonance. Thus while some new 
physics does arise at the scale of tree level unitarity violation, it is new physics described 
by the existing low energy EFT. 

Of course this is certainly not a universal cure, and if a pole arises on the physical 
sheet it can potentially lead to the emergence of a negative norm (ghost) state or failure of 
analyticity. The point is that the violation of tree level unitarity is not synonymous 
with the breakdown of the EFT or the need for new fundamental d.o.f. as is 
often assumed®. Neither does it have to be the occurrence of a single resonance that cures 
the problem, but it could be composite resonant/bound multi-particle states. 

In summary then, we can never truly diagnose the breakdown of unitary using per¬ 
turbation theory, we can only diagnose the breakdown of perturbation theory itself. If a 
theory has a Hermitian Hamiltonian, there is no logical impediment to the theory being 
unitary non-perturbatively despite the fact that it may violate tree-level unitarity at some 
scale. Some new physics does arise at the scale of tree-level unitarity violation, for instance 

®We note in passing that precisely such a resummation was proposed as a means to deal with the pertur¬ 
bative expansion of nonrenormalizable field theories within the Efimov-Fradkin/superpropagator approach 
for which the real part of the perturbative amplitude grows exponentially [36] . 

^Although in many cases it is correct - e.g. the Higgs mechanism has now been confirmed as the correct 
solution for the unitarity violation of the Proca-Yang-Mills version of the standard model. 
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a physical resonance, be it sharp or a composite metastable multi-particle state. However 
the description of this new physics may only require the existing EFT. If something like this 
happens for Galileon theories, then this resolves our apparent problem with incorporating 
the Vainshtein mechanism quantum mechanically. 


2.4 Absence of UV fixed point 

The idea that a theory can make sense of itself non-perturbatively despite being pertur- 
batively non-renormalizable is part of the asymptotic safety program [38]^®. In essence, 
one assumes the existence of a UV fixed point around which we can define the theory. For 
instance naively one may argue that for Galileon theories, at high energies we can suppose 
that the highest order Galileon operator in Sa dominates. If one such operator dominates 
then the action is scale invariant, however not under the scaling associated with the engi¬ 
neering dimensions. 

In d dimensions, if an n-th order Galileon operator dominates at high energies and this 
defines the theory near a UV fixed point, then the conformal weight S defined via 

{7r{x)7r{y)) ~ , (2-20) 

will be S = (d -|- 2 — 2n)ln. However unitarity, expressed through the need to have a 
well-defined Kallen-Lehmann spectral representation, i.e. a unitary representation of the 
conformal group, requires for scalars S > d/2 — 1, i.e. the conformal weight must always 
be greater than that for a free scalar field in the appropriate dimension [43, 44]. This 
inequality becomes 

2>n, (2.21) 

which cannot be satisfied except for the trivial case n = 2 which is just the usual free field. 
In fact in the case of the (d-|- l)’th Galileon n = d-|- 1 we shall perform an explicit compu¬ 
tation of the Wightman function in Sec. 6 and we shall find an entirely different unitarity 
behaviour with no evidence of a UV fixed point. This is equivalent to the statement that 
the Wightman functions will not be tempered distributions, i.e. their Fourier transform 
will not be polynomially bounded. 

We thus see that there is a fundamental conflict between the desire for a unitary theory 
and the existence of a UV fixed point. Our perspective will be quite different, if Galileons 
can be made sense of non-perturbatively without the introduction of additional degrees of 
freedom (by which we mean fundamental - resonance or bound state degrees of freedom are 
acceptable), then we still anticipate the length scale A to control the high energy behaviour, 
i.e. there will be no UV fixed point and the Wightman functions will not be tempered. 
This is more typical of truly gravitational theories where, for instance, both the string 
scale and Planck scale determine the properties of the high energy theory and there is no 

^°For discussions related to the present case see [39-41] and for different approaches [42] 
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sense in which these are UV fixed points. Analogous systems which are argued to be UV 
finite are Little String Theories [16] where despite being ‘field theories’ their UV behaviour 
retains remnants of their stringy behaviour and there is no UV fixed point. In fact, as 
we discuss later, our perspective is that Little String Theories are the most closely related 
to our conjectured behaviour of the UV properties of Galileons. More generally this idea 
of UV completion without a UV fixed point is referred to as ‘asymptotic fragility’ in the 
scenario constructed in [20, 21], to distinguish it from asymptotic safety. 

2.5 Classicalization, Vainshtein and Galileons 

The ‘classicalization’ proposal [23-26, 45] is another such proposal of how a theory may 
self-unit arize above its strong coupling scale A. Specihcally it proposes that the scattering 
of particles at energies S S> A is dominated by the emergence of semi-classical configura¬ 
tions, ‘classicalons’, which ultimately decay into many soft quanta [46]. The typical number 
of such quanta N scales as V ~ Er^{E), where r^:{E) is the Vainshtein/classicalization 
radius defined below. The distinction between this and the self-healing mechanism [37] is 
that the new states which emerge are not single resonances but quasi-classical bound multi¬ 
particle states. These are nevertheless not new degrees of freedom. This picture is modeled 
after the ‘known’ behaviour of gravity whereby high energy scattering will produce ‘messy’ 
black holes which in turn go through the process of balding and ultimately evaporation. 
Galileon field theories satisfy precisely the conditions argued to lead to such a phenomena. 
That such an analogous behaviour can occur for a field theory such as Galileons is less 
surprising once we view the Galileons as the helicity-zero graviton in massive gravity. We 
shall go further and argue that any theory that classicalizes must have an inherent grav¬ 
itational non-locality which is ultimately responsible for its different high energy behaviour. 

We will find in subsequent sections arguments that support this picture using the 
Galileon duality. Specifically we will find that the high energy behaviour of Galileon fields 
is dominated by their semi-classical contributions to the path integral as manifested in 
the exponentially growing density of states ~ ~ gSr*(E) £qj. Galileon operators. 

This behaviour is consistent with the emergence of configurations with many soft-quanta 
whose ‘entropy’ scales as Er^{E) where r^{E) is the Vainshtein/classicalization radius. 
Concretely the peak of the resonant part of the spectral density is from finite number of 
particle states for which N ~ Er^{E). These properties rest on a nontrivial UV/IR mixing 
due to the Vainshtein mechanism which connects high energy scattering amplitudes with 
large semi-classical configurations. We attribute these properties to a fundamental gravi¬ 
tational non-locality that arises in these theories. 

In Galileon theories, in common with the classicalization arguments [23-26], associated 
with an incoming state of energy E = ^/s we may define a radius r* which characterizes the 
region of strong coupling for scattering. This is none other than the Vainshtein radius for 
the system. To see this we note that for Galileons in four dimensions, the Vainshtein radius 
is usually given in terms of the mass of some external source which creates the background 


- 13 - 


configuration as 


ry = A 


-1 


M 

M^iJ 


( 2 . 22 ) 


However if we compute the classical self-energy in the Galileon configuration associated 
with this source we find a hnite contribution of the form 


~ r 


-1 

V 


M 

Mpi 


(2.23) 


That this self-energy is hnite is the hrst clue to the better UV properties of these theories. 
Rearranging these formula we hnd 


rv{E) ~ A ^ 



(2.24) 


which in the language of [23], is the classicalization radius r^{E) = rv{E). In a generic 
Vainshtein/classicalizing theory we may parameterize this as r^{E) = A~^ (^) with 
a > 0 (here we use the notation a to denote the order of the spectral density to be dehned 
later). In general dimensions 2a — 1 = l/{d + l) for Galileons. We shall continue to refer to 
this as the Vainshtein radius understanding that it is being expressed in terms of the energy 
of the Galileon held and not the energy of any potential source for the Galileon. It will 
become apparent that this is a more appropriate parameterization and it is independent 
of how we choose to couple to the Galileon to other helds. 


We know that classically the Vainshtein mechanism works in the sense that whenever 
an energy E is localized in a region less than r^{E) the associated Galileon held will be 
highly non-linear, i.e. strongly coupled. In the context of massive gravity, the Vainshtein 
radius is for the helicity-zero mode of the graviton what the Schwarzschild radius is for 
the helicity-2 mode, namely the scale at which nonlinearities become important. This also 
means that in massive gravity, the effective size of the gravitational held around a massive 
source, e.g. a black hole, is set not by the Schwarzschild radius, but by the Vainshtein 
radius, at least from the point of view of scattering Galileons or particles coupled to the 
Galileon held. Particles which are not directly coupled to the Galileon will continue to see 
the Schwarzchild radius as the effective size. 


Although the Vainshtein radius is usually determined for spherically symmetric sources, 
it is easy to make a scaling argument that in any system, for which an energy E is localized 
in a region of size ry, then ry is the radius which determines the onset of strong coupling. 
Note that it is the localization of energy that is important for this argument. For instance 
for planar symmetric solutions of energy E nothing special happens when scattering at 
these energy scales [47]. However, this is a consequence of the fact that in these solutions 
the energy is not localized in all spatial directions. 

In particular, this means that in scattering two high energy Galileons with incoming 
energy y/s, when the impact parameter b becomes comparable to r*(s) = ry(s) we expect 
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strong quantum non-linearities to kick in [23-26]. A similar argument is made in the 
context of GR, when two particles are scattering at an impact parameter less than or equal 
to the Schwarzschild radius associated with the incoming state r* = Rsis) then we expect 
a Black Hole to form. In the case of gravity, if the incoming particles have ultra-planckian 
energies ^/s Mpi, then the scattering will become nonlinear at a length scale which is 
much bigger than the Planck length Rs ^ The larger the energy, the larger the 

distance at which scattering is nonlinear. However this implies a UV/IR mixing we would 
not expect to arise in any local field theory. In the case of GR this is attributed to the 
fact that it is a gravitational theory, and gravitational theories have a necessarily weaker 
notion of locality. Even away from the Black Hole production region, it is well known 
that ultra-planckian scattering in the eikonel regime is well described semi-classically using 
Einstein gravity regardless of the precise UV completion of gravity above the Planck scale 
[48]. 


2.6 Loop counting parameter 

The emergence of semi-classical behaviour at high energies can be anticipated by a simple 
scaling argument. Consider the leading order Galileon interactions in d dimensions (a = 
d/2 -|- 1) 

/ di 

d'^xA^vr^a^eenV", (2-25) 

n=0 

and suppose we construct the S-matrix element between two coherent states 




J Rttc 


iS'A 


(2.26) 


where S'jj^ = Sa+ boundary terms associated with the coherent states. Associated with the 
incoming and outgoing coherent states we can define the energy by 

Ea = j (^2 + (V7r)2) |a), (2.27) 

Ef, = I (^2 + (Vvr)^) 1/3) . (2.28) 

Let us suppose for simplicity that ~ ~ E. Next we define dimensionless space-time 

and field parameters, 

x = r*x, t = rR , tt = VZtt . (2.29) 

The energy, E, dimensionally scales as 

E~y d‘^~^x y^Tr)^ • • • = (r*)'^“^Z J d‘^~^x ^971)^..., (2.30) 

and so we may identify E Z. All of the different terms in the Galileon action 

will become of the same order when ddn ~ A'^. Reexpressing in terms of dimensionless 
variables this occurs when 

(2.31) 
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Solving for Z [Z = rf A^'^) and substituting in the previous equation we find that the scale 
at which all interactions become comparable in fixed energy scattering is 

\ ( E \ ‘*+1 

r.m = X (a j ’ 

for any Galileon model. As we have already emphasized this scale is precisely the Vainhstein 
radius expressed in terms of the energy of the Galileon field rather than the mass of the 
source which might create it. This is a more appropriate representation for the purposes 
of scattering processes. 

Finally, rewriting the action in terms of dimensionless variables we have 

S'j, = EtJ'a, 

d+2 

= (^-j 5 :,. (2.33) 

S'j^ is the dimensionless Galileon action with boundary terms (i.e. effectively with A = 1 
and E = 1). This shows us that the proper loop counting parameter, i.e. the quantity 
playing the role of h is 

d+2 

Lie = (A/^)^ , (2.34) 

which decreases as energy increases. In other words, in terms of dimensionless field and 
space-time variables, the path integral is 

/ d -\-‘2 ^ 

. (2.35) 

This suggests that quantum mechanics becomes less and less relevant as energy increases 
or more precisely that the semi-classical approximation to the path integral becomes better 
and better at high energies [23-26] (for similar observations from a different perspective 
see [32]). This result is a manifestation of the UV/IR mixing implied by the Vainshtein 
mechanism. As it stands this argument is of course too simplistic since we have not 
demonstrated that however we resolve the UV divergences of the LEEFT will be consistent 
with this scaling. For instance, if we take the theory with a defined hard cutoff at Ac then 
the rescaled theory will have a hard dimensionless cutoff at r*(E)Ac. On the other hand 
perturbative unitarity continues to break down when the dimensionless energy is unity. If 
however the theory self-unitarizes by some mechanism, then this scaling will be correct. 

2.7 Including irrelevant operators, 5 b 

The previous argument was made including only those special Galileon terms in Sa and 
one may expect it do be destroyed as soon as we include the infinite number of terms in 
Sb which will generically be present in matching to a specific UV completion. However at 
this point the Galileon symmetry comes in useful! As we have already remarked, all the 
terms in 5 b have at least one additional derivative relative to those in 5^ 

5b ~ ^5a . (2.36) 
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Naively one might think that then S'b ~ {E/A)Sa and so at high energies 5b will dominate. 
However, in reality, in the semi-classical regime we should scale 


d ~ 


n{E) 


d 


(2.37) 


at least in determining the contributions to the semi-classical action which will determine 
the leading contribution to the phase shifts and absorption coefficients. In other words 

5r 1 

asi7/A^oo. (2.38) 

Thus not only do we expect that the high energy scattering to be semi-classical, we expect 
that the leading contributions to the scattering amplitudes to be determined entirely by 
the special Galileon operators Sa- This is closely analogous to the well-known behaviour 
of quantum gravity, that in certain limits, the leading contribution to the trans-Planckian 
scattering amplitude is determined entirely by the Einstein-Hilbert contribution to the 
action, regardless of the details of the UV completion [48-50]. 


2.8 Unitarization at High Energies and Non-polynomial boundedness 

How then do we expect the scattering amplitude to be consistent with unitarity at high 
energies given that it violates tree level unitarity? Consider the 2 —>■ 2 scattering amplitude 
for a generic Galileon theory. The generic form of such a scattering amplitude that is 
consistent with unitarity is 

I - °° 1 

A{s,t{cOSe)) = J _ ^'2 39 ^ 

^ ^ ^ 1=0 * 

Here Si{s) is the scattering phase shift and is the absorption coefficient reflecting 

the inelastic parts of the scattering amplitude. Unitarity is maintained provided that 6i{s) 
is real and (3i{s) > 0. 


If the conjecture is true, that at high energies the scattering is dominated by semi- 
classically configurations in the path integral, then the phase shift and absorption coeffi¬ 
cients are simply determined by the real and imaginary parts of the classical action eval¬ 
uated on classical conhgurations sourced by the scattered particles with the appropriate 
c.o.m energy and impact parameter b = l/p. 

7?e[5ciassicai] = 2(5i(s), (2.40) 

Im[5classical] = 2/3;(s). (2.41) 

The classical configurations considered are in general complex solutions of the equations 
of motion which may tunnel through classically forbidden regions. It is these latter con¬ 
tributions that will give rise to a nonzero Im[5] as in standard instanton calculations. For 
example, in the case of scattering in gravity, this imaginary part arises because even in 
the absence of other interactions the two incoming particles could form a black hole. The 
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likelihood for this to occur is determined by the Entropy for a black hole of energy ^/s and 
angular momentum I, i.e. in this case we expect 

ImfS'classical] = 2A(s) = S'entropy (-E, 0 • (2-42) 

This is the ‘Black Hole ansatz’ of [50] (see [51] for an updated review and [49] for earlier 
arguments). The remaining part of the ansatz is that Si{s) takes the form 

6i{s) = TTk{E, l)SentropyiE, l)/2 , (2.43) 

where k{E, 1) is assumed to have a milder E dependence. 


We thus see that unitarity is guaranteed at high energies if the scattering becomes 
quasi-classical provided only that Im[S'ciassicai] > 0- This restriction is equivalent to re¬ 
quiring that the Euclidean action is positive definite. Since we are in the quasi-classical 
region we may expect to be able to replace the sum over I by an integral, and then this 
integral may be performed by a saddle point/stationary phase approximation. In other 
words there will be some generically complex Z* which minimizes the combined phase shift 
and exponential. We expect this either to be at the scale for which the impact parameter 
b ~ \l*\/y/s is comparable to the Vainshtein radius, i.e. when |f*| ~ y/sr^{s) which becomes 
increasingly larger for high energies. For generic angles we utilize the large I expansion for 
the Legendre polynomials 


Pi{cos9) ~ — 


Then the amplitude will be 


V27rl sin 6 


^i(l+l/2)d+iTT/4: _ ^—i{l+l/2)d—in/4: 


(2.44) 


H(s, f(cos 0)) PS — / 

Jo 


dl 


Vl 


0 \/ 27r sin 6 


g2i(5,-2/3;+i(Z+l/2)0-|-i7r/4 _ g2j(5;-2ft-j(Z+l/2)0-i7r/4 


(2.45) 

where primes denote derivatives with respect to L In this case there are two saddle points 
for I which will depend on the angles according to 


dl + ' dl 


(2.46) 


Assuming one or the other of the two phases contributions we get in the saddle-point 
approximation. 


A(s, t(cos 6)) PS =F 


VT* 


I 


^2i6i^-2l3i^±i{U + l/2)B±i-K/A 


\/2sin0^_^5'' 


(2.47) 


Although explicitly constructing the classical conhgurations is extremely difficult, simple 
scaling arguments for Galileons in four dimensions predict the large s behaviour 


<5/*(s) ~ \/sn(s)co(6i) = A ®/^s^/®co(6») 


(2.48) 










and similarly for f3i^{s). To see this it is sufficient to consider a semi-classical configuration 
of energy ^/s and compute the action integrated over the light crossing time for the strong 
coupling region r*(s). 

Thus, rather than blowing up, we may expect the scattering amplitude for fixed angle 
9 to oscillate or decay exponentially at large s as 

A{s, t(cos 9)) ~ . (2.49) 

For generic Vainshtein/classicalizing theories we would be led to 

^(s, t{cos 9)) ~ ~ g-c(0).“/A2“ ^ (2.50) 

which is consistent since for Galileons in four dimensions r*(s) = ry(s) ~ 

[23-26]. The limit we have just described is the standard quasi-classical limit in which the 
classical scattering cross section is recovered and so we may expect the total cross section 
to grow as the square of the impact parameter which we identify with r*(s) 

/2 

o-totai(s) ~ vr— ~ vrr^(s). (2.51) 

s 

In the above, the growth parameter a in r*(s) ~ (y^/A)^“ ^ will be correlated with 

the growth properties of the Wightman functions. 

The emergence of the generically non-integer power for s, and hence A in the expo¬ 
nent, at first sight seems peculiar from the perspective of perturbation theory, but it can 
arise from the growth of certain entire functions which can be obtained by resumming the 
perturbative expansion. We shall see this behaviour explicitly in the computation of the 
Kallen-Lehmann spectral densities. This is the characteristic behaviour of a non-localizable 
field theory whose momentum space Wightman functions for Heisenberg picture operators 
grow exponentially as W{k) ~ , 

This conjectured high energy behavior can be entirely consistent with unitarity, but 
what it does violate is polynomial boundedness. The power a = 3/5 ensures that there 
will be some direction in the complex s plane for fixed angles in which the amplitude grows 
exponentially. This follows simply from the analytic properties of functions in the complex 
plane. More precisely, the conjectured behaviour H(s,t(cos0)) ~ violates the 

Cerulus-Martin bound [52] whose derivation assumes polynomial boundedness. Specifically 
the Cerulus-Martin bound states that the fixed angle scattering amplitude cannot fall off 
faster than 

H(s,t(cos 9))\>C . (2.52) 

A simple minded argument for this bound comes from noting that if the amplitude falls of 
as then on continuing s into the complex plane s —?• we find 

which for a > 1/2 will inevitable grow exponentially in the upper half complex s plane 
thus violating polynomial boundedness. Again this is not in contradiction with unitarity 
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since this occurs for unphysical momenta, but it does imply an inherent non-locality. 

The semi-classical ansatz for the high energy cross-section grows faster than allowed 
by the usual Froissart bound. For massless Galileons, or for Galileons whose intermediate 
scattering states include massless particles, this is acceptable since the Froissart bound 
does not apply. For the massive case, i.e. for which there is a mass gap for all states, this 
behaviour would be forbidden for a local field theory, but is consistent with the bounds 
for non-localizable theories which we discuss in more detail in [14]. The form (2.52) is al¬ 
lowed for non-localizable theories where the bounds on the growth of the amplitude in the 
whole complex plane are much weaker |A(s,t)| < [53] and so the Gerulus-Martin 

bound does not apply. These amplitudes may then violate the Froissart bound [53] which 
makes the explicit use of polynomial boundedness. It is important to note though that 
while this behavior violates polynomial boundedness at fixed angles, if a < 1 we never¬ 
theless expect that the fixed momentum transfer t < 0 amplitude, including in particular 
the forward scattering limit, to be polynomial bounded. This follows from the fact that 
Im(74(s,t)) < Im(A(s,0)) for real t < 0 and use of the Phragmen-Lindelof theorem. We 
will discuss these issues in more detail elsewhere [14]. 

Thus if the classicalization proposal for unitarization at high energies is correct also 
for Galleons with a mass (regardless of how small), then these theories are necessarily 
non-localizable. The non-localizability appears to be consistent with the arguments of [23- 
26, 45] where the exponential decay of the elastic scattering amplitude is interpreted in 
terms of the product of intermediate quasi-classical multiple particle states, ‘classicalons’, 
of size r*(s). As in the Black Hole example, the production of such states from scattering 
implies an inherent non-locality over the scale r*(s). Even in the massless case we expect 
this non-locality to show up in the violation of polynomial boundedness, at least away from 
forward scattering and fixed physical momentum transfer limits [14]. 

Thus, we see that there is, in principle, no obstruction to UV completing the theory in 
a Lorentz invariant, unitary and even analytic way. However we must give up polynomial 
boundedness [14]. This is the signature of an inherent non-locality which we regard as 
signalling that these theories are quintessentially gravitational as we discuss in the next 
section. 

3 Gravity and Locality 

The fundamental distinction between gravitational theories and ordinary gauge theories 
is the following theorem: ‘There are no local diffeomorphism invariant observables in 
gravity’. The theorem is simple to prove since even a scalar field transforms under a dif¬ 
feomorphism 6^(p{x) = but its implications to attempts to quantize gravity 

are profound, as was noted early on by Wigner and Salecker [54, 55]. In particular, it is 
no longer clear how to apply the notion of locality and causality in the strong quantum 
gravity regime since the usual definition of the vanishing of the commutator of fields 4>{x) 
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outside the lightcone is neither gauge invariant nor independent of the background metric 
used to define the field. The normal approach to deal with this is to focus on a gauge 
invariant quantity such as the S-matrix, and hnd a substitute for locality and causality, 
namely polynomial boundedness and analyticity. It has been argued, however, that gravi¬ 
ton scattering amplitudes do not satisfy polynomial boundedness, at least at fixed angles 
[50, 51], and that any UV completion of gravity includes an inherent non-locality [56-58] 
that could even potentially explain aspects of the information loss paradox [59, 60]. We 
briefly review the arguments in favor of an inherent gravitational non-locality 


3.1 Giddings-Lippert locality bound and UV/IR mixing 

The existence of a locality bound in gravity is also connected with the existence and 
properties of black holes. In scattering two high energy particles at high energies, it is 
expected that when the impact parameter becomes less than the Schwarzschild radius 
associated with the centre of mass energy, 

b<Rs{Vs), (3.1) 


a black hole will form. As the energy increases the black hole increases in size, indicating 
a UV/IR mixing. This UV/IR mixing in itself runs counter to the expectations of a local 
field theory for which UV physics decouples from IR and is indicative of some, albeit mild, 
degree of non-locality. In particular, Giddings and Lippert have proposed a fundamental 
bound on non-locality in gravitational theories [61, 62] related to this. To state this bound, 
consider a smeared field (l)x,p = / fx,p{x')(j){x') which effectively creates a wavepacket 

at position x with momentum p. For example this may be achieved with a Gaussian wave 
packet 

4 ,, = I . (3,2) 

In a local field theory there is no problem constructing a dressed two-particle state of the 
form \ijj) = 4>xi,pi4'x2,P2\^) where jO) is the gravitationally dressed vacuum. Furthermore, 
by the Wightman axioms, if the test functions fxx,px{x) and fx 2 ,P 2 ix) have compact support 
(which is allowed if the correlation functions are tempered distributions) and the support 
of one region lies entirely at spacelike separations for the support of the other then we 
would have for the equal time commutation relation 

[^Xi,piJ 4 'x2,P2] = 0 - ( 3 - 3 ) 


If there is overlap as in the case of the Gaussian wave packets then we would expect a 
weaker statement to the effect of 


< 


1 


,-\x-x'\y{2L^) 


nYXuPl,YX2,P2il - 

However in a gravitational theory, once the following condition is violated 


(3.4) 


xi - X2I > Rs(|pi +P2I), 


(3.5) 
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gravity becomes sufficiently strong that a black hole forms, or at least that the fluctua¬ 
tions in the geometry are so strong that it is no longer possible to talk about locality in 
the usual sense. Concretely this would imply a violation of the locality of the Wightman 
axioms through a violation of the condition (3.4). 

In practice what does this violation of locality mean? From the perspective of the S- 
matrix for particle scattering, which by the LSZ construction is related to the time ordered 
Wightman functions, we expect a violation of locality to imply a violation of polynomial 
boundedness since the latter is derived on the basis of locality [15, 63-65]. This violation 
may only show up in some unphysical directions in the complexified momentum space since 
unitarity typically constrains the growth in the physical directions. This is indeed observed 
in the eikonel behaviour of scattering amplitudes, although it is conjectured not to arise in 
the forward scattering limit t = 0 and fixed momenta transfer t < 0 [50, 51]. 

From the perspective of the Wightman functions, we will argue that the failure of 
locality shows up in the failure of one of the central assumptions of the Wightman ax¬ 
ioms, namely that the correlation functions are tempered distributions which is a sufficient 
condition to imply polynomial boundedness of the S-matrix [15, 63-65]. All tempered dis¬ 
tributions have a Fourier transform. Our conjecture for the properties of the Wightman 
functions implies that they exist in momentum space, but do not admit a Fourier transform 
back to real space, because of the locality bound. More precisely the Wightman functions 
will be argued to be GelTand-Shilov distributions. 

3.2 Quantum Mechanics of M-theory 

A closely related argument which emphasizes a different aspect of the role of black holes was 
given by Aharony and Banks [66]. Consider quantum gravity (e.g. M-theory?) with a fixed 
asymptotic form for the geometries, e.g. Minkowski or AdS. There will exist some asymp¬ 
totic symmetry group which we assume to contain at least time-translations. Quantum 
mechanics then guarantees the existence of a generator for translations, i.e. a Hamiltonian, 
H which can be used to dehne Heisenberg operators 

d{t) = e*^*6e-*^*. (3.6) 

We have been intentionally agnostic about what the operator is or even its space dependence 
to be as general as possible. As in field theory, it is natural to define Wightman functions 
and their associated Kallen-Lehmann spectral representation 

poo 

W(t) = (0|0t(t)0(0)|0) = / dHe-'®Vo(^), (3.7) 

Jo 

where the spectral density is manifestly positive 

po{E) = Y,SiE-En)md\n)\\ (3.8) 

n 

\n) are a complete set of energy eigenstates (including bound states) and the only assump¬ 
tion is that all these states have positive energy (technically bounded below is sufficient 
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since a constant part of H drops out of 0{t).) 

In order for the Wightman function to exist it is necessary that the spectral density for 
the chosen operator po{E) to grow no faster than an exponential. A typical operator will 
have a spectral density which scales with the density of states. For a field theory whose UV 
behavior is determined by a UV fixed point we expect the density of states in d dimensions 
in finite volume V to behave as p(E) ~ ^j^-g note that for gas of 

massless particles at temperature T the energy in d spacetime dimensions will be U ~ VT'^ 
and so the entropy will scale as ^entropy ~ EjT ~ Ej{E, hence the density of states 
is 

p{E) ~ ~ ^ 

This growth is sufficiently slow, i.e. slower than a linear exponential, in all dimensions 
that the Wightman function exist. There will also exist special operators whose growth 
is at most polynomial. These are the operators usually used to build the theory and its 
interactions. 


By contrast for a gravitational theory, we expect the density of states at high energies 
to be determined by the density of states for black holes which scales as 

p{E) ~ entropy = ^c{E/ 

In other words, as in the previous S-matrix argument, we expect the physics at high energies 
to be dominated by the production of black holes which are known to have exponentially 
large density of states. Contrary to the field theory case, this grows faster than a linear 
exponential in any dimension d > 3. This stronger growth is related to the negative spe¬ 
cific heat capacity of Black Holes in asymptotically Minkowski spacetime. Thus, operators 
whose spectral density scales as the density of states po{E) ~ p{E) can no longer give rise 
to well defined Wightman functions, i.e. they are no longer associated with well defined 
localized Heisenberg fields. 


The resolution is to look at operators for which the high energy behavior is suffi¬ 
ciently cutoff, meaning in practice that the spectral densities are cutoff above the Planck 
scale. This can easily be achieved by working with modified operators whose inner prod¬ 
uct between the vacuum and energy eigenstates decays exponentially at high energies, e.g. 
O' = e~^ ^ ^ however this resulting operator will not be local. Inevitably this 

implies that there is no precise notion of time locality, i.e. no tempered distribution W{t) 
for the original operators 0{t). Alternatively this high energy cutoff may be achieved by 
smearing the original operator over a length scale L comparable to the Planck scale . 


d{t) ^ d{t) 




1 

2L^ 


h-i')' 


0{t') 


(3.11) 


The Wightman function of the smeared operator is well defined 

roo 2 

W{t) = {0\d\t)d{t)\0) = dEe-^^\-'^^po{E), (3.12) 

Jo 
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provided that > 2cM^^. There is no contradiction with the validity of the EFT at 
low energies since as long as we simultaneously make observations at energies below Mpi 
and distances above the Wightman functions can be replaced by their smeared values 
which satisfy the usual requirements of a local field theory. 


3.3 Giddings-Lippert bound for time non-locality 

At a mathematical level, the claim of [66] is that the Wightman functions are no longer 
tempered distributions, however they do exist as distributions but must be defined using 
test functions from an appropriate GelTand-Shilov space [67] Sa, which contain no functions 
of compact support. If f{t) and g{t) are drawn from such a space then 

IT(/,5) = (0|Ot(/)O(ff)|0), (3.13) 

exists and is finite. The fact that we must use test functions which do not have compact 
support to define the Wightman functions already indicates that the operator is non-local. 


We may use this conjectured behaviour of the Wightman function spectral densities 
to derive an obvious generalization of the Giddings-Lippert bound. Suppose that 0{x) is 
now a space-time dependent field. Again we express the vacuum expectation value using 
the Kallen-Lehmann spectral representation 

roo p 

{0\6^{x)d{x')\0) = dgpoig) (3.14) 

Jo —_^2 

where 

^dk = J-^^e{k^)2Tr6{k'^ + p). (3.15) 

The conjecture of [66] similarly implies po{k) ~ gc(/i/Mpi) large p for d = 4. Now 
Giddings and Lippert [62] consider the equal time commutation relation, and smear these 
operators with wave packets which give a spread in space. For instance, defining the 
smeared operator 

4, = (3.16) 


then the vacuum expectation value of two equal time product of two such operators is given 
by 


(0|6t^p4'y|0) = / dppoip) , 

Jo 7/ c 2 —_^2 

= / ^0(fe°)d(-A;2)27rpo(-A:')e'"-("-"')e-^-^. (3.17) 


It is transparent that despite the spatial smearing the k^ integral does not converge. Thus 
the non-locality implied by the pure spatial form of the Giddings-Lippert bound applied at 
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equal times is not sufficient. We must extend the bound to include an inherent non-locality 
in time 

\ti - t 2 \ > Rs{\pi +P2\ + \Pl +P 2 I) • (3-18) 

Of course in a Lorentz invariant theory it is not surprising that non-locality in space implies 
non-locality in time, nevertheless we see that the finiteness of the Wightman functions has 
more to do with smearing them in time. To demonstrate this consider now the time smeared 
operators 

dx,E= [ dr ^ . (3.19) 

J—00 y/Z tt Jj 

Then we have 


e^x\E'\^) = / 

' Jo 


' k?‘ = — u? 




Now convergence of all the integrals is guaranteed by the factor e = e . 

Reorganizing we have 




(e'^+e'^)lP f°° 


I d/j, po(p)e / dke 

0 J/c2 = —^2 




In order for the fi integral to converge we need L? > c/Mpj. We see from the form of the 
integral that there is a spread in the time dependence of the order 


A(t - t') r.{E + E')E >{E + E'YjMl^ > Rs{E + E'). 
This is the time non-locality version of the Giddings-Lippert bound. 


(3.21) 


3.4 Effective Field Theory point of view 

These argument are sufficiently simple that they are meant to apply to any UV completion 
of gravity in all dimensions greater than 3. As we discuss in the next section, a field whose 
spectral density grows slower than a linear exponential are known as strictly localizable 
fields, and those whose spectral density grows faster than a linear exponential are known 
as non-localizable fields. Hence the claim is that one of the key distinctions between a local 
held theory and a gravitational theory is: 


In a gravitational theory, the spectral densities of generic operators grow 
faster than linear exponentials, i.e. operator valned fields in quantnm gravity 
are fnndamentally non-localizable. 


All of this is consistent with the original observation that no gauge invariant local 
observables exist in gravitational theories. How can we see this from the point of view of 
the EFT description of gravity? Consider again quantum gravity, treated now in an EFT 
sense, for spacetimes whose asymptotic symmetry group includes time translations. As long 
as we are at low energies we may express the metric as g^uix) = g^v{x) -|- where 

g is the background solution with the desired asymptotics. Then follow some gauge hxing 
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formalism for which the constraints may be solved locally in time, for instance the BSSNOK 
formalism. By construction this formalism is designed so that the solution of the constraints 
only requires inverting elliptic spatial operators. This reduces the system to the physical 
degrees of freedom in a way which remains local in time. We may then solve the Heisenberg 
equations of motion for the remaining physical degrees of freedom perturbatively using the 
Yang-Feldman formalism [68] or by one of the many other approaches. The Yang-Feldman 
approach has the virtue for example of making causality and locality manifest at each 
order in perturbation theory. Even if solving the constraints involves inverting the 
Yang-Feldman equations for the propagating degrees of freedom always take the schematic 
form 

Ahp{x) = ( 3 - 22 ) 

where A is some hyperbolic operator of the schematic form H ~ □ and Tp{x) is the 
effective stress tensor source including non-linearities in hp. This may then be solved in a 
local/causal way using the retarded propagator 

hp{x) = h%{x) - j <1%Gr{x - y)^fp{y). (3.23) 

where Gp is the retarded propagator for A, and Up{x) is the quantized free physical field . 

We can now compute the spectral densities of some time-local operators (but not 
necessarily space-local) 0{t) = F{hp{x)) perturbatively as an expansion in E/Mp\. At 
any finite order in perturbations there are a finite number of counterterms that need to be 
added to remove divergence, nevertheless matching against some specific UV completion 
will determine specific values for these counterterms. The spectral density will now be 
expressed as a perturbative expansion (with possible logarithmic terms from loops) 

Since the EFT is local, and the operators defined are local in time, at any finite order 
in perturbations the Wightman functions will exist as tempered distributions. Similarly 
at any finite order in perturbations the S-matrices will be polynomially bounded. In this 
language passing to the UV completion simply corresponds to summing the series to an 
infinite order. If this sum converges absolutely and Pq{E) grows slower than a linear 
exponential then the Wightman function would remain well dehned. But this would imply 
that it was possible to define gauge invariant observables which are local in time (even 
if they were non-local in space due to the solution of the constraint equations) which we 
know is impossible from time diffeomorphism invariance alone. Even if this series were 
asymptotic then assuming it can be made sense of by Borel summation or some similar 
procedure, as long as it grows slower than a linear exponential we would be left with the 
same contradiction. Thus the absence of local gauge invariant observables can only be 
consistent if Pq{E) grows as fast as or faster than a linear exponential. This may be easily 
achieved even if the coefficients c„ are on average positive and fall of at a particular rate. 
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For instance if Cn ~ jnl in d = 4 we would recover the prediction Pq{E) ~ ec(_E/Mpi)2^ 
Of course we cannot from an EFT point of view determine the c„ precisely because of 
the need to add counterterms, they can only be determined by matching against a specific 
UV completion, nevertheless we can at least see how the locality of the EFT of gravity 
can be consistent with the non-locality of the UV completion based on the convergence 
properties of the perturbative expansion. The effective field theorist that truncates the 
expansion at small finite order will see no evidence of non-locality, but the intrepid field 
theorist that computes to large orders, and can perform the EFT matching, will begin to 
see the non-locality by the convergence properties. 

3.5 The case of Little String Theory 

Above we have reviewed suggestive arguments that imply that in gravitational theories, the 
spectral densities of certain gauge invariant observables grow exponentially fast, consistent 
with the requirement that such operators cannot be local. Operators associated with these 
observables are known as non-localizable fields and we shall review the definition of these 
in Sec. 4. Another class of theories in which this behaviour has been argued to occur are 
Little String Theories (LST’s) [69, 70]. LSTs describe M-theory compactihed on a and 
are thus six dimensional theories. They are obtained as decoupling limits of N coincident 
fivebranes in string theory in which Mpianck —)• oo, 0 keeping the string scale Mgtring 

fixed. Since Mpianck —^ co they are no longer gravitational in the usual sense, and so one 
might be tempted to think that they are field theories. Indeed, the infrared limit of LSTs 
are normal six dimensional superconformal field theories (SCFTs), but the full LST cannot 
be a local field theory (i.e. there is no UV fixed point) because it includes non-trivial string 
effects since Mstring is kept finite. As a result we expect this theory to maintain certain non 
field theory features of string theory - hence the name ‘Little String Theory’. In particular 
it is expected that it preserves T-duality. This alone implies that LSTs cannot be local 
field theories since the T-duality transformation is itself non-local. In [16] it was proposed 
that LSTs are quasi-local theories, i.e. non-localizable field theories of the Jaffe type with 
a = 1/2 (see Sec. 4 for the definition). 

The salient feature which establishes this interpretation is that the Wightman functions 
of LSTs can be shown to have the exponential growth characterstic of the Hagedorn density 
of states [66, 71, 72]. In particular the Kallen-Lehmann spectral density for the two-point 
function is expected to grow as 

p{p) ~ e (3.25) 

for an order unity constant c, where N is the number of five-branes. A similar growth 
is expected to arise for all correlators, and this growth cannot be removed with a simple 
wavefunction renormalization. Our perspective will be that LSTs are the closest known 
analogue theories to Galileons in the sense they are local field theories in the infrared limit, 
but they are fundamentally non-localizable in the UV, having no UV fixed point, and an 
exponentially growing density of states. Both LSTs and Galileons are obtained as double 
scaling limits of gravitationally theories. Both theories exhibit a duality symmetry which 
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is a remnant of a symmetry of the original theory they are defined as decoupling limits of. 
This analogy is summarised in the table below. 


Theory 

Little String Theory 

Galileon 

Parent Gravitational Theory 

M-theory compactified on T® 

Massive Gravity/DGP etc. 

Decoupling Limit Theory 

LST 

UV Galileon 

Infrared Limit 

Six Dimensional SCFT 

Galileon LEEFT 

Decoupling/Scaling Limit 

9s ^ 0, A/pianck ^ ^ 

m 0, Mpianck OO 

Strong Coupling Scale 

A = Mstring/V^ 


Infrared Region 

E <C Mstring/\/iV 

E <C A 

Symmetry of Parent Theory 

M-theory T-duality 

Diffeomorphism Invariance 

Duality of D.L. Theory 

T-duality 

Galileon Duality 

Growth of Spectral Density 

lnp{ij.) ~ cMstving^/p/VN 

G GG „(d+2)/(2(d+l)) 

lnp{p) ~ ^(^+2)/(d+i) 

‘Vainshtein’ Radius 

n{E) = \/]VM,"kg 


Jaffe Class 

Quasi-local a = 1/2 

Non-localizable a > 1/2 


Since LSTs are argued to be quasi-local, the analogue of the Vainhstein/classicalization 
radius is independent of energy scale. In [16] it is argued that LSTs should be understood 
as Jaffe type quasi-local field theories and we shall follow an analogous description in the 
case of Galileons in what follows. 

3.6 Galileon Duality and the Locality Bound 

We have previously reviewed the arguments to suggest that gravitational theories, and 
certain decoupling limits of M-theory/string theory, exhibit some degree of ‘gravitational 
non-locality’. This non-locality will lead to a violation of the usual Wightman axioms, 
specifically the temperedness assumption. This in turn leads to a violation of the assump¬ 
tions used to derive the polynomial boundedness of the S-matrix, about which we shall 
discuss more elsewhere [14]. 

In the following we conjecture that these same properties hold for Galileon theories, in 
other words Galileon theories are intrinsically gravitational. Our main concrete evidence 
for this proposal is that we are able to give an explicit off-shell UV complete quantization of 
a specific Galileon model, and demonstrate that its Wightman functions grow exponentially 
as 

p{E) ~ (3.26) 

where r^{E) = A“^(£'/A)<^+1 is the Vainshtein/classicalization radius. This exponential 
growth means that there do not exist local off-shell observables for Galileons, although 
there do exist precisely defined non-localizable off-shell fields. This exponential growth 
means that we can only talk about locality for separations satisfying the analogue of the 
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Giddings-Lippert locality bound 


1 ( E \ <^+1 

|x|»r,(ii;) = -(^-j . (3.27) 

Mathematically, this corresponds to stating that Galileon fields are no longer ‘operator 
valued tempered distributions’ but rather ‘operator valued Gehfand-Shilov distributions’ 

[67]. 


Although we only demonstrate this explicitly for one specific example of a Galileon 
model whose S-matrix is trivial, we conjecture that the same qualitative exponential growth 
of the off-shell Wightman functions is true for all interacting Galileon models. This is rea¬ 
sonable since all Galileon models have the same scaling of r^:(E) regardless of the choice of 
coefficients, and the phenomenology of the Vainshtein mechanism is qualitatively the same 
in all such models. Furthermore, a semi-classical argument reproduces the same conjec¬ 
tured behaviour for all Galileon models regardless of choice of the coefficients. 

This exponential growth of the correlation functions will feed into an exponential 
growth of the scattering matrix in some directions of the complex plane (albeit constrained 
by the requirements of analyticity and unitarity). Thus, we conjecture that the interacting 
Galileon models necessarily violate polynomial boundedness of their scattering amplitude 
[14]. 


Another viewpoint on the origin of the locality bound is provided by the existence 
of the duality between various Galileon models which we review in Sec. 5. At the level 
of the field theory Lagrangian, the duality is an equivalence map, enacted by a non-local 
field redefinition, between two naively distinct Galileon theories. However, as we review in 
Sec. 5 from the point of view of massive gravity where the duality was derived in [73, 74] 
(see also [75]), or from the coset construction considered in [76, 77], the duality is simply 
a choice of diffeomorphism gauge. In other words, by viewing the Galileon models as a 
gauge fixed version of a diffeomorphism invariant theory, which is the natural interpreta¬ 
tion both in the massive gravity and coset perspectives (see Sec. 5), the Galileon duality 
becomes a symmetry, just as T-duality of LST is a symmetry of M-theory. Assuming that 
the quantization respects this symmetry, which seems desirable, then following the usual 
arguments we would led to conclude that there are no local off-shell gauge invariant op¬ 
erators in Galileon theories. Of course we may choose to fix unitary gauge to go back to 
the original Galileon Lagrangian, but then we would be left with the statement that there 
are no local off-shell operators invariant under the Galileon duality transformation, which 
is the remnant of the previous statement. This is analogous to the statement that in LST, 
there are no local off-shell operators invariant under T-duality, which as we have argued is 
one of the clues that LSTs are not local field theories. 

Thus in our proposed picture, the Galileon duality transformation, far from being 
an accident, is viewed as a central property of the Galileon theories which enforces the 
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locality bound. Furthermore the duality must hold at the quantum level, since it is simply 
a remnant of a local symmetry of the UV completion (assuming that the latter is not 
anomalous). We shall in fact show that the same exponential growth of the Wightman 
functions occurs in two different duality frames, albeit for two different operators, consistent 
with the duality map. 

4 Strictly Localizable versus Non-localizable field theories 

The interest in field theories with a weaker notion of locality than standard perturbatively 
renormalizable field theories dates back to the 1960’s, largely to two programs. On the 
one hand there were attempts within the axiomatic approach to quantum field theory to 
determine the minimal requirement for incorporating locality into the axioms of field theory. 
A precise condition was given by Jaffe [78, 79] and Meimam [80] that distinguished (in 
Jaffe’s language) so-called strictly localizable field theories, and non-localizable field 
theories. On the other hand there was an interest, sparked by Guttinger [81-83], Efimov 
[84-86], Fradkin [87], Volkov [88-90], and others [36, 91-97] to give resummation techniques 
to field theories that were non-renormalizable from the standard point of view, but could 
potentially be partially resummed into a renormalizable field theory. This approach was 
sometimes referred to as the ‘super-propagator’ [98-103] or ‘Efimov-Fradkin’ [104-106] 
method and at that time the most plausible application was to non-renormalizable theories 
of the Weak nuclear force. The resummed operators generically fall into the class of fields 
with an exponentially growing spectral density. It was even argued that this approach may 
be applicable to gravity [107-109] and potentially tame the naive non-renormalizability of 
gravity. This original program failed and was ultimately replaced by the now standard field 
theory tools of renormalization group and effective field theories. In particular the modern 
perspective is that all non-renormalizable non-gravitational field theories should be viewed 
as EFTs. It will be our contention here however, that Galileon theories are intrinsically 
gravitational and are explicit realizations of such non-localizable field theories. 

4.1 Definition of Strictly Localizable, Quasi-local and Non-Localizable fields 

The physical meaning of the definition of a strictly localizable versus a non-localizable field 
is easy to understand by considering its associated Wightman functions written in the form 
of a Kallen-Lehmann spectral representation. Gonsider some operator 0{x), its position 
space Wightman function may be defined as 

poo 

W{x,y) = {0\d\x)6{y)\0) = dy,po{fi)Wf,{x,y), (4.1) 

Jo 

where poifJ') is the Kallen-Lehmann spectral density. Unitarity is encoded in the statement 
that po{l^) > 0 for all /r > 0 and is real. Stability is encoded in the statement that poifJ-) 
only has support for p, > 0. The function W^{x,y) is the Wightman function of a free 
scalar held of mass 
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A strictly localizable field is one for which the spectral density integral f dfi con¬ 
verges so that W{x,y) can be given meaning for x ^ y without smearing (i.e. using test 
functions of compact support). In the case of four dimensions W^{x, y) has the asymptotic 
properties 

W^{x,y)p:i (4.2) 

W^{x, y) ^ -^e-W4 ^ _ y)2 < 0 , ^|x - yp » 1 ,(4.3) 

(47r|x — y\y' 

and a similar exponential behavior occurs in general dimensions. It is clear that the con¬ 
vergence of the spectral density integral requires that the spectral density grows no faster 
than a linear exponential in yfji, i.e. that 


or equivalently that 


lim y/Ji-^\upo{lj) = 0 
At->-oo dy 


' Inpoiy) 

dy - 7rn ^— < oo , 


/ 4 ' 


3/2 


(4.4) 


(4.5) 


There are obvious generalizations of this condition to the n-point Wightman functions. If 
W{{ki}) denotes the Fourier transform of the n-point Wightman function as a function of 
complex momenta then we require the bound 


W{{ki}) < 1*^*1 


(4.6) 


Whenever these conditions are satisfied for all n-point Wightman functions, the operator 
is said to be strictly localizable. 


A non-localizable field is then any field for which 


lim Jy-^ Inpo{y) = oo 

11^00 dy 

(4.7) 

/■“ lnpo(/^) 

1 ^3/2 

(4.8) 


When this condition is true the position space Wightman functions do not exist, i.e. we 
may not use test functions of compact support. However, the Wightman functions do 
exist as we shall see when defined using an alternative class of test functions, namely those 
from the appropriate GelTand-Shilov space [67]. In this case the Wightman functions are 
still distributions but not tempered ones. The description of such distributions is well 
developed [67]. A non-rigorous but practical way to understand this is to note that the 
momentum space Wightman functions are well defined, and so provided we choose test 
functions which are sufficiently delocalized in position space, then their Fourier transforms 
will fall off sufficiently rapidly at large k to compensate the exponential growth of the 
spectral density. Explicitly we can define the smeared Wightman function by 

r°° r d'^h 

W{f,g) = J^ dy J j^^f*{k)g{k)2Trpo{-k‘^). (4.9) 
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Motivated by the theory of entire functions, we define the ‘order’ of the spectral density o. 
by the requirement that 



lim 

At^oo In y/Jl 

(4.10) 

and the ‘type’ a by 

lim 

(4.11) 



In other words, the high p behaviour of the spectral density is taken to be 



Po{lA ~ X subdominant terms . 

(4.12) 


It is then clear that a suitable choice of momentum space test functions are the set of 
infinitely differentiable functions for which f{k),g{k) < as \k\ —)• oo for some 

finite constant C. This gives a (non-rigorous) definition of the appropriate Gehfand-Shilov 
space [67] which is usually denoted as Sa- When the test functions are drawn from such a 
space, the smeared Wightman function is well-defined. The Fourier transform of such test 
functions define the dual space S'®. Their properties are well understood [67] and, roughly 
speaking, correspond to functions which are no more localized than a function of the form 

(4.13) 

Test functions of compact support are thus excluded. An obvious simple example is the 
case a = 1 for which the set of momentum space test functions fall of as Gaussians and 
hence by the usual Hardy uncertainty principle argument that position space functions 
cannot be localized more than a Gaussian. 

This discussion should be contrasted with the ‘traditional’ case, tempered localiz- 
able helds, which are operator valued tempered distributions for which the spectral density 
grows at most as a polynomial. For tempered fields, only a finite number of subtractions 
need to be performed to give the position space Feynman propagators meaning. In the case 
of non-localizable fields an infinite number of subtractions, i.e. an infinite number of renor¬ 
malizations should be performed. This connects with the notion that non-localizable field 
theories have to do with perturbatively non-renormalizable field theories whereas localiz- 
able helds have to do with renormalizable held theories [91]. Precisely one of the historical 
reasons for interest in non-localizable helds was to give a potential renormalizable descrip¬ 
tion of perturbatively non-renormalizable held theories. In particular we note that the 
usual notion of operator product expansions is only applicable in the case of tempered lo- 
calizable helds for which the short distance behaviour of a product of two, operators A[x) 
and B[y) diverges by no more than an inverse power oi \x — y\. This behaviour is only 
applicable if {ilj\A{x)B{y)\il}) is a tempered distribution in (x — y) and if an analogue were 
to exist in the case of non-localizable held it would have to be substantially different [93]. 

Finally, as a special case we can consider a quasi-local held for which 

lim 

/i^OO 


\apo{y) = hnite ^ 0 . 


(4.14) 
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Figure 1. Schematic form of the spectral density for a theory with a mass gap in each of the 
possible cases. In addition, bound states may exist between the pole and the beginning of the 
continuum. In the massless limit, applicable for the examples computed in the text, the continuum 
begins at the position of the pole located at p. = 0. 


is a limit of a non-localizable field with a —)• 1/2. This case is clearly special since the 
position space Wightman function does exist at sufficiently large space-time separations. 
As a consequence, such fields retain some degree of locality even though they violate strict 
micro-locality/micro-causality. Interestingly, it is precisely this case that is expected to 
be applicable to Little String Theories [16] where the exponential growth of the spectral 
density is associated with the Hagedorn behaviour for a string. 


In summary the nature of the field is detemined by the order a of the spectral density. If 
a > 1/2 the field is non-localizable, a = 1/2 it is quasi-local, 0 < a < 1/2 strictly localizable 
and if a = 0 it is a tempered localizable field (see Fig. 1). Textbook renormalizable field 
theories correspond exclusively to the latter case a = 0. 


4.2 Functions of a free field 

As an illustrative example of a strictly localizable field whose Wightman function is not 
tempered, consider a free field operator (/(x) which for simplicity of calculation we take 
to be a massless field. Define the operator 0{x) to be the normal ordered :: exponential 

[100, no. 111] 


CXD Yi 

d{x) = - 1 : . 

n=i 


(4.15) 


- 33 - 









An elementary application of Wick’s theorem tells us that the two-point Wightman function 
for 0(x) is given by 


{0|d(l)d(!,)|0) = ^ A^{0|:<A(ir;:.j(!,r:|0> 

nlmi 

n,m=l 

n=l 

= g9^(O|0(a:)<A(2y)|O> _ ^ _ 


(4.16) 

(4.17) 

(4.18) 


For a massless field in four dimensions the Wightman function is given by 

w{x, y) = mmym = ,- j —gi ■ ('‘■w) 

The Wightman function is a ‘generalized function’ or ‘distribution’ [67] so even products of 
it must be interpreted with great care. The problem with the eq. (4.18) is that in effect it 
includes the exponential of a delta function from the support of the Wightman function on 
the light cone, which is meaningless as it stands. To deal with this more rigorously we must 
smooth the operators with test functions to give meaning to the operator product. However, 
to proceed it is simplest to directly compute the spectral density via the perturbative 
expansion. Concretely we have 

271 /* 1 

e{k^)27rpo{-k^) = y (4.20) 

OQ 2n 

= e{k^)2^5{-k^) + y ye{k^) VLn{-k^) (4.21) 

n\ 

n=2 


where Qn{—k‘^) is the n-particle phase space density 


n^i-k^) = 

This is given by (see Appendix A) 


~ n ^ 

P[ J dki 

. 1=1 


{2y6^^\k -yki). 


i=l 


1 




n—2 


(4.22) 


(4.23) 


(Idvr^)*^ ^(n — 2)!(n — 1)! 

The perturbative summation converges to give a pole contribution plus an entire function 


27rpo{li) 


(167r2)"-n!(n — 2)!(n 
2vr5(l^) + (; 2, 3, 


which behaves asymptotically as /i —t- 00 as 


Poikl ~ 6 24/3,r2/3 


1 

7! 


(2yy^ 


.n—2 


1)! 




(4.24) 

(4.25) 


(4.26) 
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We thus see that 0{x) corresponds to a strictly localizable field with 0 = 1/3 and has 
faster than polynomial growth. Unitarity is preserved since p(/r) > 0 for all /i as is evident 
since each coefficient in the entire series is positive. Note that since the continuum contri¬ 
bution to the spectral density is an entire function, there is no ambiguity with regards to 
its non-perturbative definition, i.e. there is no need for Borel summation. 


An example with a stronger high energy growth is given by the operator 

(j){x) 


0{x) =: 


1 - g(t>{x) 

The position space Wightman function is naively 




\n-l-l . 


(4.27) 


{0|6(i)6(s)|0) = 

n=0 


+ 1 )! 

(4^2)n+l ((f_y )2 


1 

((xO - 2/0) - ie)2)«+i ■ 


(4.28) 


This series is asymptotic and so ill-defined without further specification. A natural guess 
would be to try to Borel resum the series, however there is no guarantee that such an 
approach preserves unitarity. 

The ‘resolution’ which ensures unitarity is to directly compute the spectral density as 
before 


2'Kpo{g) 


2tt6{h) -b ^ -b 1)! ^n+iig ), 

71=1 


2vr5(M) + £ 5^” (I67r2)^(n - 1)! 

n=l 

, , jV -b 327r^) 




(4.29) 

(4.30) 

(4.31) 


Unlike the position space correlation function, the spectral density is convergent (it is in 
effect a double Borel transform of the position space correlator - see Appendix B) for 
all g. In other words this series gives an unambiguous non-perturbative definition of the 
Wightman function. The asymptotic form of the spectral density is 


2ttpo{p) 


g 167r2 


4 

9 9 
2567r4 ’ 


(4.32) 


hence this is an example of a non-localizable field with a = 1. We see that the failure of 
the position space expansion to converge is now consistent with the fact that the field was 
not strictly localizable. One may be concerned that 0{x) is ill-defined due to the lack of 
convergence at 4>{x) = Ijg, however we can resolve this by making g complex and taking 
the limit to real g. Since the spectral density is entire, it is automatically the boundary 
value of its analytic continuation. There is no non-perturbative ambiguity in determining 
the spectral density. 
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Remarkably, this exponential behaviour is possible to anticipate based on very simple 
semi-classical considerations which will anticipate what we find using the Galileon duality. 
The strong coupling region for the operator 0{x) is the region ~ I/5. Consider a 
semi-classical field configuration which is strongly coupled inside a region of radius r*, i.e. 
the Vainshtein/classicalization radius. The energy of the field in this region will scale as 

(4.33) 

J\x\<rt 9 

Turning this around, the Vainshtein/classicalization radius for this operator is 

n{E) ^ Eg'^ ^ . (4.34) 

Finite contributions to the action (e.g. from complex instantons) will then scale as ~ 
Er^{E) ~ fj.g'^. The spectral density thus grows as 

p{g) ~ ~ ^ . (4.35) 

This example nicely anticipates our general theme, that the Jaffe order of growth, a, is 
determined semi-classically. 


4.3 Higher n-point functions 

The construction of higher n-point functions proceeds analogously by using the appropri¬ 
ate spectral representation form for the n-point function [112-114]. For example, for the 
exponential function 0{x) = fn" • — 1 : we have the 3-point function 


00 00 00 ^ 2 {n-\-m,-\-p) 


(0|d(xi)O(x2)O(x3)|0) = EEE 




where 


1 


n=l m=l p=l 

1 


n!m!p! '•12-13-23: 


47r2 [xi — XjY — {{x^ — x^) — fe)2 


We then express this in Fourier space as 

/■ j4 


where the 3-point function spectral density is 

9^ 


p{kl2,ki3,k23) = EEE 


n=l m=l p=l 


nlmlpl 


'^n{~ki2)^m{~ki^)i^p{ — k2^) 


(4.36) 


(4.37) 


(0|6(xi)0(x2)0(x3)|0) = /^/^/(4.38) 


(4.39) 


and it is understood that the integration is over the range < 0 and k^j > 0. 

As a consequence of the exponential form of the interaction we have in this case 


p(/ci2, fel3) ^ 23 ) — {‘^'^)^Po{~ki2)Po{~ki3)po{ — k23), 


(4.40) 
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where po{—k‘^) is the 2-point function spectral density (4.24). We thus see the spectral 
densities are still entire functions and that the high energy behaviour of the 3-point function 
in each —kfj is the same as that for the 2-point function. This implies 

(0|0(xi)0(x2)0(x3)|0) = (0|6(xi)0(x2)|0)(0|6(xi)0(x3)|0)(0|6(x2)0(x3)|0) . (4.41) 

which was to be expected given the exponential form of 0(x), however we should remember 
that this expression only has meaning after smearing with appropriate test functions. 

4.4 Canonical Quantization of Vainshtein theories and non-localizablilty 

Using the duality map, we will argue below that Galileon theories correspond to non- 
localizable quantum field theories. It is tempting to conjecture that this property is true of 
all theories exhibiting the Vainshtein/classicalization mechanism. This connection is based 
on the fact that in these theories, the Vainshtein/classicalization radius r^{E) grows with 
energy which means that the (Euclidean) action on associated semi-classical configurations 
grows with energy faster than a linear growth S ~ Er^(E). This simultaneously means 
that the semi-classical approximation is better at high energies and then we may anticipate 
that the spectral density of a typical operator is determined semi-classically as p ~ e'^. The 
order of growth of the spectral density is then 

1 1 d In r* 1 

ry = — -/> — 

2 2dlnE 2’ 

implying that the associated operators are non-localizable. Turning this around, any UV 
embedding of a theory which exhibits the Vainshtein mechanism for which the UV theory 
is a local field theory will inevitably destroy the Vainshtein mechanism since locality would 
enforce r* to decrease at high energies so that p grows slower than an a linear exponential. 
This is consistent with known attempts [115]. 



We will now see that the non-localizability of Vainshtein type theories can be antic¬ 
ipated from a naive application of the rules of canonical quantization in the interacting 
representation. The S-matrix may be formally defined in terms of the interaction Hamil¬ 
tonian density l-Lmtix) as 

5 ^ ^ (4.43) 

For the S-matrix to be Lorentz invariant, we require that this notion of time ordering is 
Lorentz invariant. This requires that 




= 0 , 


{x-yf >Q. 


(4.44) 


In order for this to be the case the field T-imtix) must be a strictly localizable held otherwise 
the product of operators in position space has no meaning. 


In the interacting representation Tdmt {x) is built out of the free helds which themselves 
do commute outside the light cone, and so a sufficient condition for Ihii-atix) to be strictly 
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localizable is that it is built out of a finite number of powers of fields and their derivatives. 
However, as we have seen we can also have strict localizability even with an infinite number 
of powers of the field, as exemplified by the exponential example 'Hintix) =: : — 1. 

The characteristic feature of theories exhibiting the Vainshtein mechanism is that their 
Hamiltonian is built out of an infinite number of powers of fields in such a way that 
the Hamiltonian density is not a strictly localizable field. The origin of this is that the 
Vainshtein mechanism requires a modification of the kinetic term. A finite number of local 
corrections to the kinetic term in the Lagrangian corresponds to an infinite number of 
terms in the Hamiltonian as a result of the Legendre transform between them. 

4.4.1 Anti-DBI 

As a simple example of this phenomena consider the DBI Lagrangian with a possible switch 
for the signs in the kinetic term 

C = eA^ — + e(5^)2/A"^. (4.45) 

The ‘right sign’ DBI theories is e = +1 and the ‘wrong sign’ anti-DBI theory e = —1. It 
is know that for e = +1 the low energy S-matrix is consistent with the possibility of a UV 
completion with a local theory [ 11 ] and in fact it is easy to find local completions using a 
two field system (e.g. [116]) or as a brane moving in higher dimensions. On the other hand 
for e = — 1 there is no local UV completion^^. 

Let us nevertheless proceed naively and construct the interacting Hamiltonian in order 
to canonically quantize this theory 

Uint = -eA^ + eA^ : 1 + 1 + e{Vct>f (V^)" :, (4.46) 

where P is the momentum density conjugate to </>. Due to the square root structure this 
Hamiltonian contains an infinite number of terms and so we should check whether T-Lint {x) is 
a localizable field or not. In fact it is non-localizable which may be easily seen by analyzing 
each square root separately. Defining 

°° prQ/ol ^n+l 

d{x) =: e^l + eP2(a;)/A4 - ^ ^ . (4 47 ) 

^ ^ ^ ^r[ 3 / 2 -n]r[n + l] A^*" ^ ^ ^ ’ 

then we have the Wightman function 

(0|6(i)6te)|0> - E (4.48) 

^ (r[ 3/2 - n]r[n + l])^ A«^ 

The spectral density is given by 

9(t»)2ipo(fc) = + f; (r[3/?/!|]r{r+l])4 A^”^"<'^> 

we construct the two field model as in [116] the second field would be a ghost and so there is no 
unitary UV completion. If we construct the brane EFT one of the extra dimensions would be time-like and 
this would lead to violations of causality and unitarity. 
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^r!(n —r)! (27r2)’^ 
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r=0 

n 

E 
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r!(n-r)! (27r2)- ^ g^ 2 (n-r-) ^)- 


(4.50) 
-(-xg)("-’-) (4.51) 
37, (4.52) 
(4.53) 
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(x 2 — (x^ — ie)2)3"--»’ 


The asymptotics of the spectral density at large k^ » |/c| are then 


^.W^E . 


Q‘ 2 (n-r) 


r=0 


r!(n — r)! (27r^)"- 


dk, 


2 {n—r) 

0 


f^3n—r(^o) , 


(4.54) 


D^{k) = ktr-^Y. 
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!(_l)r4n-r(4p)3n-r(g^ - 2r - 4)! 


^ r!(n — r)!(27r2)"^(167r2)3^-»’-i(3n — r — 2)!(3n — r — l)!(4n — 4)! 


,4„_4 2’^7r3/2-2«(3n-5/2)! 


!?.(/.) = C-^ ( 4 n- 4 )!( 3 n-l)! ^ ^^^ 

At large n this may be approximated as 


L4n —4 


Dn{k) = akQ 

where a ~ 0.85. With this approximation we have 
2TTpo{k) = 2'K5{k^) + 


2\kQa 
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9A;8 


81927r2A8 
which at large /cq behaves as 


(4.55) 


(4.56) 


(4.57) 


6’3’8’2’3’6’ ’ 167772167r4A8 J 


9 128\/27r 34/3 

2'xpo{k) ss 2n5{k^) H - -4 —e8,r2/3 a'I/3 , 


(4.58) 


This exponential growth is precisely what is expected from semi-classical arguments. The 
Vainshtein radius for this system is r* = A“^ [E /and so the semi-classical action 
scales as S ~ . An analogous calculation may be performed for the full interac¬ 

tion Hamiltonian with a similar result. 


From an EFT point of view we need not be concerned about this behaviour since we 
should treat the interaction perturbatively where at each finite order of perturbations T^int 
remains local. More precisely, in an EFT field theory we can add to the original Lagrangian 
an infinite number of increasingly irrelevant operators consistent with the symmetry which 
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will contribute to increasing powers of fe/A in the spectral density computation. Thus we 
cannot take this exponential growth of the leading order operators too seriously since it 
is possible that the additional irrelevent operators cancel the high energy growth. In the 
case of the DBI model, e = +1 we expect this to be the case since we know we can UV 
complete this model in a local way, e.g. via a two-field system or by viewing it as the 
LEEFT of a brane moduli in extra dimensions. For the anti-DBI case e = — 1 there will be 
no local UV completion, and this suggests that in this case there is no way to resolve the 
non-localizability of the interaction, and that this means that this theory is fundamentally 
non-localizable. However, it is also known that it is only in the case of e = — 1 that there 
is an active Vainshtein mechanism when this field is coupled to sources. This suggests that 
the Vainshtein mechanism and the non-localizability go hand in hand. For different reasons 
similar conclusions were reached in [34] with regards to the classicalization proposal. 

4.4.2 Goldstone model 

A similar model which is simpler at the level of the Lagrangian but more complicated at 
the level of the Hamiltonian is the Goldstone model where the Lagrangian is just 

£ =-1(50)2+ ^(d<^)4. (4.59) 

As in the DBI case the Vainshtein radius for this system is r* = A“^ (F/A)^^^ and so we 
anticipate that the spectral density scales as ^ . Unlike the DBI case there are only 

a finite number of terms in the Lagrangian, nevertheless there are an infinite number of 
terms in the interaction Hamiltonian which takes the inelegant form 

H = 1 : P 2 : +i : : ( 3^,2 + _ y^f) ,, ( 4 . 6 O) 

where v is the solution of the cubic equation 

P = v + j^y{v^-{V(f>f). (4.61) 

As before, we may compute the spectral density of this interaction and we find, due to the 
cube root structure of the solution for v, the predicted exponential growth in the spectral 
density. Once more in the ‘right sign’ case, e = -|-1, we would regard this as harmless and 
resolvable by adding irrelevant operators to the Lagrangian. In the wrong sign e = — 1 
case, which is the case which exhibits the Vainshtein mechanism, we expect to be stuck 
with the non-localizability. 

4.4.3 Cubic Galileon 

Let us now consider how this argument applies to the relevant case of the cubic Galileon 
model. Here the Lagrangian is very simple 

^ = -^(5vr)2(l + ^°vr) ■ (4.62) 
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and the Hamiltonian interaction takes the form 


^int = : 


2 ■ 1 + V^tt/AS 


: + 


: {V7r)^V^Tr : . 


2A3 


(4.63) 


The second term is manifestly local and harmless. The ‘dangerous’ part of the Hamiltonian 
comes from the operator 0{x) = ■ Following the same procedure if we compute 

the spectral density of this operator we find that it grows exponentially as po ~ ^ 

which is the expected behavior given that the Vainshtein radius is r* = A“^ [E/. 


Unlike in the previous two examples there is no sign issue since the sign of the inter¬ 
action is degenerate with the field redefinition vr —?■ —tt. Both signs are bad in that it is 
known that the cubic Galileon does not admit a local Lorentz invariant UV completion 
[11]. Consequently, and consistent with the picture developed from the duality described 
below, we anticipate that the cubic, and indeed all Galileon models, are non-localizable 
field theories. 


One interesting difference with the Galileon model is that even if we attempt to view it 
as an EFT, it is known that because of the Galileon symmetry all the additional operators 
5 b that we can write down that are Galileon invariant have additional derivatives. This 
means that the contributions to the spectral density from these additional operators always 
have a different form than those from the leading order operators. In the semi-classical 
regime the additional suppression amounts to a suppression by 


d 1 _ My/® 

A "" M 


(4.64) 


If we assume that the semi-classical approximation to the spectral density is good at high 
energies then we find that it will be dominated by its prediction from the leading operators 
Sa with negligible corrections from Sb- In this sense, the Galileon symmetry prevents any 
possible additional irrelevant operators from resolving the non-localizability of the leading 
order ones which is consistent with the fact that the Galileon theory does not admit a 
local UV completion. This is of course not a rigorous argument because it could be that 
the semi-classical approximation is not good (despite naive expectations based on the path 
integral scaling), nevertheless it sits well with the known properties of Galileons. 


4.4.4 Implications of a non-localizable Hamiltonian 

If, as we have argued, it is the case that all Vainshtein type theories have a non-localizable 
interaction Hamiltonian then we are led to the following consequences: 


• There can be no quantization of a Vainshtein type theory for which the fields are 
localizable. If there were it should be equivalent to the canonical procedure which 
requires only locality in time. 

• We can maintain the canonical quantization, i.e. causal and unitary evolution in time 
and micro-causality, at the price of breaking Lorentz invariance. 
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• If we insist on maintaining Lorentz invariance, the quantization will not be canonical 
non-perturbatively (although at any finite order in perturbations it will be canonical) 
since there will be no local evolution in time. 

To reiterate the last point in a different language, and to set ourselves up for our future 
perspective: we can potentially quantize Vainshtein type theories, e.g. Galileon models, 
in a Lorentz invariant unitary way, but we cannot maintain micro-causality (at least with 
its usual definition) since that would require a well-defined operator that generates time 
translations which we have excluded (we cannot maintain Lorentz invariance without that 
operator being localizable). At least for one Galileon model, albeit one with a trivial 
S-matrix, we shall be able to provide precisely such a quantization procedure using the 
Galileon duality transformation described in the next section. 

5 Classical Duality Map 

In this section we review the essential details of the Galileon duality discovered in [73-75] 
(see also [117]). We will make use of this duality in the next sections to provide a non- 
perturbative definition the Galileon theory which is dual to a free field. However, here we 
shall only be concerned with the classical properties of the map. 

The origin and role of the duality can be viewed in a number of equivalent com¬ 
plementary ways, as a field dependent diffeomeorphism, a field redefinition, a Legendre 
transformation, a change of representation of the Galileon coset, or equivalently as a choice 
of diffeomorphism gauge in the coset construction. We briefly review these different ideas, 
and also discuss how the duality extends to coupling to matter as was first observed in 
[ 118 ]. 

5.1 The Duality as a coordinate transformation 

The Galileon duality transformations are a one parameter family of invertible field redef¬ 
initions. Given a field 7r(x) we can define the dual field 7r(x) via the implicit relations^^ 
[74] 

f —)■ x^ = x^ -|- ^d^7r(x), 

Vs -.{ _ ^ . (5.1) 

[ (p^{x) = d^7r{x) —> iPfiix) = df,Tr{x) = (p^{x) 

Here A is a fixed energy scale and s is the parameter of the group transformation and 
(T = d/2 + 1 where d is the number of spacetime dimensions. Below we shall give explicit 
formula that are more practical to use than these implicit relations for actually determining 
the dual fields. The implicit relations are however practical to determine the transformation 
properties of the action. The inverse of the transformation is = P_s, 

. (5.2) 

[ d^7r(x) —> 5^7r(x) = (9^7r(x) 

choose a slightly different sign convention as in [74] so that the fields are equivalent when s = 0, 
7r(x) = 'k{x) + 0{s). 


- 42 - 



These implicit relations can equivalently be written as 

Vs : 7r(x) —^ it(x) = 7r(x) + ^{dT:{x)f , 

V_s ■■ ^{x) —> 7r(x) = fr(x) - ^{d^{x)f . 

The infinitesimal form of the duality transformation is 

67r{x) = n{x) - 7r{x) = -^{dTT{x)f + 0{s ^). 
We see that the infinitesimal transformation is a local held redehnition. 


(5.3) 

(5.4) 

(5.5) 


5.2 Duality Group 

The duality map forms a continuous group which can be extended to the full group GL{2, R) 
considered in [77]. We shall however focus on the only truly non-trivial part of the group 
determined by the above transformation laws. Performing a second duality transformation 
with parameter s' starting from n{x) and denoting the new dual held as 7r(x) we have 


x^ 


Vs 


+ ^5^7r(x) 


7r(x) 7r(x) = 7r(x) + ^ ^57r(x)j = 7r(x) + {dTT{x) f 

where we used the relations (5.1) and (5.3). This leads to a combined transformation 


(5.6) 


Vs' oVs = Vg+s' ■ 


(5.7) 


In other words the duality map forms an abelian group with transformation law s” = s + s'. 
Again we note that this group transformation leaves invariant the derivatives of the Galileon 
helds 


dfj,TT{x) = df,n{x) = df,n{x). 


(5.8) 


Dehning the Galileon invariant combination 


1 


n/^,(x) = —<“a„5,7r(x), 


(5.9) 


and similarly I{^u{x) = d^diy^{x )/and with the hat variables then we have in matrix 
notation the relations 


n-i = u-^ + si, 
fi-i = n-i + s' I, 
n-i = n-i + (s + s') I. 


(5.10) 

(5.11) 

(5.12) 


A particularly useful result is that 

det[l + sn]det[l — sH] = 1. 


(5.13) 
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5.3 Covariant Coset construction 


Here we review the coset construction for the Galileon duality [76, 77] emphasizing the 
duality as a choice of diffeomorphism gauge. Our starting point is the Galileon group 
Gal(d,l) which has been shown to be the natural starting point for the coset construction 
of the Galileon models [119]. The full set of generators are momenta Pa, Lorentz trans¬ 
formations Mab, Galileon boosts Ka and Galileon translation (shift symmetry) C. The 
Galileon algebra is (we use physics convention which includes the factor of i ) 

[Mab, Mcd] = i (vacMbd - VadMbc - VbcMad + VbdMac) (5-14) 
[Mab, Pc] = irjAcPB - ir]BcPA (5.15) 

[Mab, Kc] = ivacKb - itibcKa (5.16) 

[Pa,Kb] = iVABC (5-17) 

[C,Pa] = [C,Ka] = 0 (5.18) 

It is transparent that in this algebra Pa and Ka play the same role. Furthermore we may 

take linear combinations of Pa and Ka to define the dual Galileon boosts Ka = Ka — sPa 

with the algebra remaining unchanged under the replacement Ka -a- Ka [76]. This ambi¬ 
guity is the origin of the duality transformation. However rather viewing it as an automor¬ 
phism of the algebra, we shall reinterpret as a change in coordinates, i.e. diffeomorphism 
gauge. This point of view was emphasized in [77] and more clearly connects with the in¬ 
terpretation of the duality in massive gravity. 


When dealing with the breaking of spacetime symmetries, in defining a coset it is 
usual to split the generators into the broken and unbroken generators giving special role to 
the unbroken translations (for a recent discussion and references see [120]). For instance 
for the current breaking pattern Gal{d, l)/ISO{d — 1,1), Galileon boosts Ka and shifts 
C are broken, and rotations Mab unbroken. Translations generated by Pa are assumed 
unbroken. The coset element is then parameterized in an asymmetric way in which we 
include the broken generators and the unbroken translations 

5 ^ ^ ^iPAX^^iPx)C^ ( 519 ) 

where 7r(x) and 4>^{x) are the Goldstone fields from the broken Galileon shifts and boosts. 
However in the specification of any space, a given point x spontaneously breaks transla¬ 
tions. Thus it is more symmetric to view that all translations are always broken, in other 
words the coset is G/H = Gal{d, l)/SO{d — 1,1) so that H is the isotropy (not isometry) 
group. This is precisely what we mean when we specify Minkowski spacetime as the coset 
Minkoswki= ISO{d — 1, 1)/ SO{d — 1,1). 


Since all translations are being viewed as broken, we express the coset element in terms 
of only broken generators with associated Goldstone fields 


^ — giPAY^ix) ^i<f>^{x)KA^iTr{x)C ^ 


(5.20) 
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In this picture Y^{x) are the Goldstone fields for the translations which are viewed as bro¬ 
ken by the need to define a point x^. In other words these are the coset analogues of the 
Stiickelberg fields in the massive gravity/bi-gravity construction of the duality [73, 74, 118]. 

Since id is a one-form, it is manifestly invariant under diffeomorphisms even though we 
have not chosen to gauge translations. This is as it should be since the coset formalism, as 
is well known, separates the concept of gauged translations from diffeomorphisms. Thus 
even though there is no ‘gravity’, we must demand that the actions we construct respect 
diffeomorphism invariance and this gives any theory of a nonlinearly realized spacetime 
symmetry a gravitational flavor. 

The Maurer-Cartan one-form is given by 

— = (IY^Pa + d(f)^KA + {4>AdY^ + d'x)C 

= ojpPA + ^o^Ka + wcG (5.21) 

where we have introduced the one-forms uip,u:^,LOc for future notational convenience. In 
this form Pa and Ka, and hence Y^ and (f)^ are treated symmetrically. Indeed defining 
^ = IT — Y^(j)A we have equivalently 

- in-^dn = dY^PA + d^^KA + i-YAdcj)^ + dTT)C (5.22) 

which is equivalent to (5.21) under the interchange Y,P —cj), —K. 

Legendre transformation 

The next step in the Galileon coset construction is to impose the inverse Higgs constraint 
[119] 

= (t’AdY^ + dn = —YAdcj)"^ + dn = 0 . (5.23) 

These equations can be formally solved as cj)"^ = —dnjdYA and Y^ = —dnjdcl)^ so that 

TV = TT — Y^d-KldYA (5-24) 

TT = t: — , (5.25) 

which are just a Legendre transform and its inverse. This confirms that the coset con¬ 
struction naturally reproduces the Legendre transform picture for the duality described in 
[75]. 

Galileon Wess-Zumino actions 

The one-forms u}p,uj^,uic are the natural building blocks from which the Galileon La- 
grangians are constructed. The resulting action should be diffeomorphism invariant and 
Lorentz invariant. Since we are imposing the inverse Higgs constraint ojc = 0 then naively 
we can only construct the Lagrangian directly from ojp and o;^. However, it transpires 
that all wedge product terms are total derivatives and all other terms will lead to higher 
order equations of motion. The resolution is to construct the Galileon Lagrangians as 
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Wess-Zumino terms [119] by introducing an auxiliary (d+l)-th dimension and defining the 
Lorentz-invariant (d+1) co-cycles 

(5.26) 

for n = 1... d -|- 1 and we have used a shorthand notation where it is understood that all 
the Lorentz indices are contracted with the Levi-Civita symbols so for example in d = 4 
^2 = ^ABCD A Up A Up A Up A u^. The An are all exact 

An = d^n , (5.27) 

where 

f^n = elT{^UKT-\^upY+^-'- - (5.28) 

2(d -|- 2 — n) 

To see this it is enough to use dup = dux = 0 and so 

dPn = edp{^UKT-\^upY+^-^ - ^!i^lI^0^d<^^6(Au;i,)’^-2(Awp)'='+2-n^ (5 29) 

and noting that after a little rearrangement [119] 

{n-l)ct>^d(t)Ae{^UKT-‘^{^upY+^-'^ = -{d + 2-n)(t^^dYAe{^UKT-^{^upY+^-'^. 

(5.30) 

whence djdn = e (tt -|- (j)'^dYA){/\UK)'^~^{f\upY~^^~'^ = An- Since the An are exact, the 
Wess-Zumino terms can be written as 

Sn= ! An= [ /3n. (5.31) 

J M‘i+^ J M‘i 

These define the terms that arise in the action Sa- 

Diffeomorphism gauge choice 

Let us reiterate that all of the above formula are manifestly diffeomorphism invariant and 
we have made no choice between Pa or Ka about which generators represent translations 
and which Galileon boosts. To write the action in standard field theory language we must 
choose a gauge. Since both Y^ and (/>^ are Lorentz vectors, it is natural to choice a 
one-parameter (s) family of gauges in the form 

Y^{x) + s (f>"^{x) = , (5.32) 

where x^ are the coordinates on the manifold. Implicit in this gauge choice is the iden¬ 
tification of Ka = Ka — sP^ as the generator of Galileon boosts since PaY^ + Ka4>^ = 
Pax^ + Ka4>^- However it is important to stress that by phrasing this as a diffeomorphism 
gauge choice, we have disassociated ourselves with the identification of who is the correct 
generator of Galileon boosts or translations. It is clear that any choice for s is equivalent 
and since this is just a gauge choice they must lead to equivalent representations of the 
physics. This is exactly how the duality is viewed in the context of massive gravity and 
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bigravity [74, 118, 121], where it amounts to a one-parameter choice of gauges for the 
Stiickelberg fields + sd^TT. Furthermore, as shown in [118], the duality easily 

extends to the coupling to any matter helds simply by writing the standard diffeomorphism 
invariant Lagrangian for matter coupled to (j)^ and and their derivatives in a diffeo¬ 
morphism invariant way. 


In this gauge the inverse Higgs constraint implies 

dn + (j)A{dx"^ — sd4>^) = 0 . 


(5.33) 


Defining if = vr — ^s4>^(f>A this amounts to dif -|- <j)Adx^ = 0 which is easily solved as 


(t)A{x) = -d^TT{x) , 


(5.34) 


if is the dual Galileon field to vr under a hnite s transformation defined through the duality 
map if = TT — ^sd^Trd^TT. We note in particular that 


f dn dY^\ f Ott d4)"^\ dir 

V5?/ ^ dy ) ^ ^ dy J dy ’ 


(5.35) 


where y is the auxiliary (d-l-l)-th dimension, so that the Galileon terms simply become (up 
to a constant) 

‘‘ ite{Au;pf-^{AuKT-K 


Sn = 


I Mi 


It thus remains only to express up and uk in terms of if 

Up = dx^ + sdd^7r{x) = {6p + sd^dB^{x))dx^ , 
uj^ = —dd^TT{x) = {—d^dB^{x))dx^ . 

Defining Up = Up + sup we have Up = dx"^ and so 
Sn = 

Sn = 


(5.36) 


(5.37) 

(5.38) 


/ pe{A{uf,-suif+^-^{AuKr-\ 

(5.39) 

JMa 


(d+l-n)! ^ 

/^ J 1 t M 1 f Jn+r ) 

^ (a -|- 1 — n — r)!r! 

r=0 

(5.40) 


where Sn denotes the usual Galileon terms with the simple replacement vr —?■ if. This is 
precisely the duality map identified in [74, 75]. This is consistent with [76, 77]. 

Let us reiterate that the value of s is purely a gauge choice. Provided there is no 
diffeomorphism anomaly, it is impossible for any physics to depend on the precise value of 
s. In this sense all representations of the coset, i.e. choices of s, are equivalent. 
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5.4 Duality in the presence of matter 

The duality transformation easily extends to the case of Galileons coupled arbitrarily to 
matter [118]. This result is clear once we recognize the duality as a diffeomorphism, be it 
in the coset language or the massive gravity/bigravity language. The matter Lagrangian 
and the coupling to matter remain local after the duality transformation. 

In this article we will be mostly concerned with two specific duality frames. For this 
reason we shall set the group parameter s = 1 from now on. The hrst frame will consist of 
a ((i+ l)’th order Galileon field tt which is dual to a free field coupled arbitrarily to a scalar 
matter field y. This frame will naturally be called the (vr, y) frame. The second frame will 
consist of the free field p (dual to tt) and its associated scalar matter field y (dual to y) 
which we will call the (p, y) frame. Below we review how to transform from one frame to 
another. 


(vr, y) Duality frame 

The Galileon action which is dual to a free field can be conveniently written as 


-^ivT, y] = / d'^x 


det(l + Ii{x)){dTi{x)f + £m(7r(x), y(x)) 


(5.41) 


where /im(7r, y) represents any arbitrary coupling the Galileon vr may have with a matter 
field y. The above action is the system we wish to quantize and we will do so by motivating 
a definition of the quantized version using the Galileon duality. 


(/?, y) Duality frame 

To transform (5.41) into the (p, y) frame we notice we use the basic duality map 

x^ = x^ + ^9^7r(x) 

Ao- 

- ^d^^p{x) 
y(x) = y(x) 

and defining = dy^dypj then using the identities 

J d'^x det(l + n(x)) = J d'^x 
J d'^x = J d'^x det (1 — S(x)) 

{dTT{x)f = {dp{x)f 

in the action we have 


(5.42) 

(5.43) 

(5.44) 

(5.45) 

(5.46) 

(5.47) 


'S'[vr,y] = J d'^x 


-- det(l + n(x))(57r(x))^ + £m(7r(x), y(x)) 


./A 

= S[p, y] . 


1 


1 1 


+det(l - S(x))£„ (^p(x) -- —(5p(x)) ,y(x) 
-^(5p(®))^ + det(l - S(x))£,n ^^(x) - ^^(5p(x))^y(x) 


(5.48) 
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In the absence of matter we simply have 


5*Galileon[7r] — j d X 


-- det(l + n(x))(57r(x))^ 


= / d“'x 




= 5fr, 


(5.49) 


6 Galileon Wightman functions 

In a local field theory, the set of position space Wightman functions encodes all the in¬ 
formation about the system. As we have argued, in theories which exhibit gravitational 
non-locality, the position space Wightman function may simply not exist in position space, 
without appropriate smearing, or at best may exist only in regions of spacetime separations 
(the quasi-local case). This is the price that must be paid for the exponentially growing 
spectral density associated with gravitational non-locality. However what may exist, and 
be calculable as we now show, are their Fourier transforms, i.e. momentum space repre¬ 
sentation of the Wightman functions. The momentum space Fourier transforms can be 
encoded in a set of spectral densities by a simple generalization of the Kallen-Lehmann 
representation for the two point function. 

In this section we will provide an explicit computation of this spectral density for the 
case of the {d -I- l)’th Galileon which is classically equivalent to a free theory. We will 
provide an operational definition of the quantum theory, and show that this dehnition is 
meaningful in a basis of a complete set of momentum eigenstates. 

6.1 Closed formula for the Duality Map 

The implicit nature of the duality map as defined in the previous section makes it unclear 
how to perform the duality transformation beyond perturbation theory. We shall now 
present a useful set of closed formulae for the duality map which can be used to give a 
non-perturbative definition. We shall make explicit use of these formula in what follows in 
defining the quantum theory. 

To begin with let us consider how to construct the field 7r(x) out of p{x). We can 
define tt{x) by its Fourier transform as 

/ Adu 

^-(fe)e*^- (6.1) 

We may then determine 7r(/c) from the inverse Fourier transform 

7r(/c) = J d'^?/e“*^'^7r(?/). (6.2) 

The trick is to now perform the duality transformation within this integral. According to 
the above map this is then simply 

n{k) = J d^y det (1 - S(y)) . (6.3) 
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Changing the dummy integration variable y to y we obtain 


vr(x) = II U{p{y)) , (6.4) 

where 

U {p{y)) = det (1 - S(y)) (^p{y) - ^ {dp{y)f^ . (6.5) 

This formula is now explicit unlike the previous duality formula. Furthermore it is non- 
perturbative, there is no need to perform an expansion in inverse powers of A. 


We can perform similar transformations to infer x from x ^md p. Specifically 

X{x) = x(x) = x(x - A~^dp{x)) (6.6) 

Switching the dummy label x —)• x we have the simpler formula 

X{x) = x{x- A-^dp{x)), (6.7) 

or in terms of the Fourier transform of x(x) 


Xix) = 


(27r)“ 


-ik.dp(x) 

e Xk e A" 


( 6 . 8 ) 


-ik.dp{x) 


We see that both of these formula contain the crucial exponential factors e a'^ which 
will be responsible for much of the quantum properties. 


6.2 Formal definition of Galileon operators (vr, x) frame 

Our aim is to define the quantum theory of a (d + l)’th Galileon which is dual to a free 
theory. At the perturbative level this map is unambiguous and a direct perturbative com¬ 
putation of the S-matrix demonstrates that the {d + l)’th Galileon theory is indeed free. 
However, within the loop diagrams we will nevertheless need to, in a given quantization 
scheme, add an infinite number of irrelevant counterterms in order to verify this statement. 
These counterterms arise from the measure of the path integral and are typically power 
law divergent. An analogous statement is that in attempting to define the Galileon opera¬ 
tors, we must choose a very specific ordering of the operators to maintain the equivalence 
with a free theory. We shall now see that normal ordering is sufficient to achieve this and 
remove all UV divergences. In this way, by the simple act of normal ordering, we may 
give a formal definition of a unitary quantum theory whieh is classieally equivalent to the 
{d + 1) ’th Galileon theory. This will be equivalent to the theory obtained perturbatively 
via canonical quantitzation with a specific choice of counterterms but non-perturbatively 
it is not canonical. 


The construction method proceeds as follows: If the {d+ l)’th Galileon theory is dual 
to a free theory at the quantum level then it must occupy the same Hilbert space, the 
Fock space of a free quantum field % = F. Given this, our goal is to define an operator 
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7r(x) which acts in this Hilbert space, and in the classical limit is equivalent to the classical 
Galileon field. Any such operator must be built out of creation and annihilation operators 
acting in the Fock space which are used to define the dual field p{x). In other words the 
operator tt{x) must be built out of the operator p{x). 


At the classical level the duality map in d dimensions can be expressed as 


{x) = j I -^Uipiy)) e' 




(6.9) 


We thus define at the quantum level the following operator 



( 6 . 10 ) 


where :: denotes the normal ordering with respect to p’s creation and annihilation operators 
and 



( 6 . 11 ) 


In choosing this normal ordering we are automatically discarding self-contractions, i.e. 
tadpole diagrams, when we come to define correlation functions. Of course, from a pertur¬ 
bative point of view this corresponds to an infinite number of subtractions. It is in this step 
that we have implicitly defined the UV theory which will contain no additional divergences, 
at least for those quantities that have physical meaning. This procedure is analogous to 
the definition of vertex operators in string theory which are also normal ordered. However, 
a fundamental difference is that the string vertex operators will have a fixed conformal 
dimension whereas these operators do not exhibit any conformal behaviour at high ener¬ 
gies. Our definition of the UV will not be associated with any UV fixed point of an RG 
flow. This is reflected in the fact that the high energy behaviour strongly depends on the 
dimensionful scale A. This is the ‘asymptotic fragility’ idea [20]. 

The operator dehnition (6.10) is further justified by the fact that if we evaluate this 
operator between two coherent states for the p field, then we will automatically get the 
classical configuration for vr which is associated with the appropriate classical configuration 
for p 


(/3|7r(x)|a) = -Kc{x) , 


( 6 . 12 ) 


where 



(6.13) 


and 



(6.14) 


In other words this definition corresponds to working in the coherent state representation 
for the path integral and performing the duality transformation of the coherent state vari¬ 
ables. 












We stress again that this definition of the operator is only formal, since we do not expect 
the position space correlators to exist without appropriate smearing. More precisely we will 
find it is well-defined in perturbation theory, which amounts to expanding the exponential 
out, but becomes ill-defined non-perturbatively without the process of smearing. The 
correlation functions are tempered distributions at finite order in perturbation theory, but 
in the non-perturbative limit they are no longer tempered. With this in mind we may 
define a regulated local field as follows^^ 


ttn{x) 



^ pik.{x-y) 
(27r)‘^ 




71=0 


(6.15) 


This operator is a well-defined local field and interactions built out of this operator will 
satisfy the conventional axioms of local field theories, including polynomial boundedness 
and temperedness of the correlation functions. Our concern then is in taking the limit 
N —>■ oo. We will see that this limit is ill-defined in position space without smearing, but 
well-defined in momentum space for the Wightman functions. 


The non-existence of the position space Wightman functions may be traced to the fact 
that the vacuum correlation of two of these fields will formally take the form 

(0|fr(rt)i(7;')|0> = / d's / d''/ j e"-'-!-) j (6.16) 

(0| : U{p{y)) :: U{p(y’)) : |0), (6.17) 

which by translation invariance of the vacuum expectation value {k' = —k) is 

(0|7r(x)7r(x')|0) = J J : [/(p(y)) :: U{m) : |0), 

= :: e=^ : |0) (6.18) 

where the function F{y, k) is a polynomial in k coming from the contractions of the ex¬ 
ponential with U and two Us together. Since this part contains only a finite number of 
terms and does not affect the convergence of the integral we ignore for now its precise form. 


Using the elementary result that 

( 0 | : :: : | 0 ) = e-^ki^^-dpd,{o\p{y)p{om 

and the expression for the d dimensional Wightman function for a free massless field 

(P(3;)P(0)) = (^_2)(|f|2_ (^0_ie)2)(-i-2)/2 ’ (6-20) 

^®This is in effect the approach proposed by Taylor to consider locality for non-localizable field theories 
[ 122 , 123 ]. 
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where 


ad = ^-^-^7T-^/^r{d/2 - 1 ) 


( 6 . 21 ) 


we find 


(0|^(x)7r(x')|0) 




where by l/yl we mean 


^ pik.(^-x'-y)rp( u\ - dtj,d^{0\p{y)pi0)\0)/ 

{27tY 

ad (Xk-yf^ 

A2. yyd_ yd+2 ^ 


d'^k 

(27r)‘^ 


^ik.(x—x' — 


y'>F{y,k) 


exp 


( 6 . 22 ) 

,(6.23) 


1 _ 1 
V- — { y ^ — ie ) 2)‘^/2 


(6.24) 


However this expression is meaningless as it stands since for any dimension d > 1 regard¬ 
less of whether k is spacelike or timelike, there is always at least one direction in which 
the integrand grows exponentially and this cannot be compensated for by the polynomial 
behaviour of F{y, k). Indeed, this expression is equally ill-defined in the Euclidean, and so 
we may not Wick rotate and use the Osterwalder-Schrader reconstruction, nor define this 
function as the boundary value in the complex time plane as usual in axiomatic formula¬ 
tions of the Wightman functions. 


By contrast, the computation of (0|7rAr(x)7rAr(x')|0) is well-defined and gives rise to a 
specific generalized function/distribution which may be extended to the Euclidean. This 
suggests that we should look for quantities that are well defined in the limit N ^ oo. It 
is conceivable that (0|7r(x)7r(x')|0) may be defined by analytic continuation, as was con¬ 
sidered often in the context of the Efimov-Fradkin/superpropagator method. While it 
is possible to analytically continue this integral into a finite integral, this mathematical 
process will in general contradict the physical requirement of unitarity. Since our aim is 
to construct a manifestly unitary theory we regard such mathematical manipulations as 
unacceptable. Indeed many of the problems associated with this historical approach to 
non-renormalizable field theories can be attributed to using analytic continuation in a way 
which breaks unitarity. 


The clue to making progress is to note that the Fourier transform of the operator 
7r(/c) = f 7r(x) is well-defined when evaluated between two multi-particle mo¬ 

mentum eigenstates. To see this we need to use the basic fact that for every Heisenberg 
operator we can write p(x) = e“*'^'^/5(0)e*^'* and similarly d^p{x) = e“*'^'^5^p(0)e*^'*. 
With this in mind we may evaluate TT{k), for which 

7r{k) = j : U{p{x))e fXX (6.25) 
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between two momentum eigenstates as follows 


{Pf\T:{k)\Pi) = j : U{pix)) : \Pi} (6.26) 

= J : C/(/5(0)) (6.27) 

= j : U{piO)) e'^ : \Pi) (6.28) 

i(Pi-Pf).dp(0) 

= {27Tf5\k - {Pi - Pf)){Pf \ : U{m) e -- : \Pi} • (6.29) 

This last expression is perfectly well-defined between any two multi-particle states. If both 
states contain a finite number of particles then the contribution from the exponential trun¬ 
cates at finite order implying that this function is polynomial bounded in k. 


In other words the following limit is well-defined 


(P/|7r(A:)|P,) = hm (P/|^jv(A:)|Pi). (6.30) 

N^oo 

Since taking the limit —>• oo is precisely the limit of resumming the perturbation series, 
then we may say that {Pf\Tt{k)\Pi) exists as a meaningful non-perturbative expression. 
Furthermore if \Pi) and \Pf) correspond to states with a finite number of particles, then 
since the expansion of the exponential truncates, {Pf\7r{k)\Pi) grows at most as a polyno¬ 
mial in k. Thus we may always perform the Fourier transform, to give a well-defined (in 
the distributional sense) expression for 

/ Adp i(Pi-Pf).ap(o) 

j^P'^-^Pf\nik)\P,) = : U{m) e -- : \P^) ■ 

(6.31) 

This is consistent with the fact that we already know a well-defined expression for (/3|7r(x)|a) 
TTcix) which may in turn be inferred from the double sum over complete sets of momentum 
eigenstates as 

{I3\n{x)\a) = '^{/3\Pn){Pn\'k{x)\Pm){Pm\a). (6.32) 

n,m 


6.3 Formal definition of Galileon operators {p, x) frame 


It is straightforward to extend the discussion to the p, x frame. At lowest order in perturba¬ 
tions in the coupling between p and x, we can treat p as a free field. Then we can compute 
the matrix elements of x from the matrix elements of y using the operator analogue of the 
duality map (6.8) 



— ik.dp(x) 

: e 


(6.33) 


which in momentum space is 


Xk 


d'^x 


^ i(k'—k).x 

(27r)^ 


Xk' 


— ik'.dp{x) 

: e 


(6.34) 
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Evaluating between two momentum eigenstates we have 


^ /* 7 /* H At -fi / 7 \ —ik^.dp{x) 

{Pf\UP^) = J : \P,] 


= / (Tx 




(27r 




-i(k+Pf-Pi).dHo) 

= {Pf\xk+{Pf-Pi) ■ e : la; 


(6.35) 

(6.36) 

(6.37) 


Once again, this expression is well-defined, and when evaluated between finite particle 
number states is polynomial in the momenta (i.e. only a finite number of terms in the 
expansion of the exponential contribute). 


6.4 Operator Products and Spectral representations 

As we have seen there is a well-defined expression for {Pf\7r{x)\Pi) , meaning that we know 
all the matrix elements of the operator 7r(x) in the momentum basis, which may be used as 
an operational dehnition of this operator. The departure from standard local field theories 
comes from trying to dehne operator products. In ordinary field theories we can make 
sense of such operator products via the operator product expansion (OPE). The product 
of two operators at points x and y may be expressed as an expansion in (x — y) in which 
there are a hnite number of leading singular terms. Divergences may be removed with a 
finite number of counterterms in the Lagrangian. In the present case since the operator 
7r(x) is built out of an infinite number of powers of p{x) this entails an infinite number of 
subtractions. Furthermore there is no conformal behaviour at high energies since there is 
no UV fixed point. These fields then do not have a standard Wilson OPE, a reflection of 
the fact that they are fundamentally non-localizable.^^ The theory never becomes confor¬ 
mal or loses knowledge of the dimensional scale A. 


In the language of axiomatic field theory, we are no longer dealing with operator 
valued tempered distributions which is an essential assumption of the standard Wightman 
axioms. To give meaning to the integral that formally defines the two-point function we 
can introduce test functions f{x) and g{x) and dehne the smeared two-point function 


(0|7r(/)7r(5r)|0) = 

// ^^'"'^Piy’k)nk)g{k)e^p 



{k ■ yf 
d+2 


(6.38) 


where 7r(/) = f d‘^xf(x)'k(x)and similarly for g. As we have discussed, the test functions 
must be drawn from the appropriate GelTand-Shilov spaces Sa- However, while the pro¬ 
cess of smearing is useful to give a well-dehned mathematical meaning to this expression, 
it disguises the physics and, in particular, the unitarity of this expression. Furthermore 
the above integrals are difficult to perform as they stand. 


note in passing that Isham, Salam and Strathdee proposed a generalization of the Wilson OPE to 
theories of the strictly localizable type in [93]. 
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A more practical approach is to recognize that since (P/|7r(x)|Pj) is well-defined, we 
may define the product of two operator by inserting a complete set of momentum eigenstates 
to define the vacuum correlation as 


(0|7r(x)7r(y)|0) = ^(0|7r(x)|Pn)(i^nK(?/)|0) . (6.39) 

n 

Again this expression is formal because the sum over intermediate states does not converge. 
However we may diagnose this lack of convergence by using the technique of spectral 
representations. Following the standard derivation of the Kallen-Lehmann (K-L) spectral 
representation we may rewrite this as 

foo !• jdu 

(0|7r(x)7r(y)|0) = y dfx J p{-k^)2TTS{k‘^ + (6.40) 

/■OO r Jd—lu 

where ojk{p) = ^Jk'^ + p and we have defined the spectral density in the usual way by 

9{k^)e{-k^)2TTp{-k^) = Y,{2T,Y5^<^\Pn - A:)|(P„|7r(0)|0)|2 . (6.42) 

n 

Unitarity requires that p{p) > 0 which will be manifest at any hnite order in perturbation 
theory. In particular, since |(Pn|'ff(0)|0)| exists and the sum over momenta is finite for a 
given pL then there is no problem computing the spectral density p{p) as we shall do below. 
Thus the content of the two-point function is contained in the well-defined finite spectral 
density p{p). This spectral density determines the Fourier transform of the correlation 
function 

(0|7r(A;)7r(A;')|0) = {27r)'^6'^{k -|- k')2'Kp{—k ^), (6.43) 

which is well-dehned. The spectral density will have an exponential growth which prevents 
us from performing the Fourier transform back to position space necessary to dehne the 
position space correlations, but we hnd no difficulties directly working in momentum space. 

Despite the absence of position space correlators, we may dehne the operators smeared 
by test functions from the appropriate GelTand-Shilov space as 

(0|7r(/)7r(5)|0) = y ^^f{k)g*{k)2Tip{-k‘^). (6.44) 

Once again this expression is well-dehned when the test functions decay fast enough at 
large \k\ to balance the exponential growth of the spectral density. This expression then 
encodes a class of well-dehned observables in this theory. 


The most general class of observables will be the set of Ai-point correlation functions. 
These may be dehned in an analogous way by dehning their generalized K-L spectral 
representations which are well-known [112-114] and then dehning the observables either 
directly in momentum space, or by using smeared functions 

(0|7r(/i)7r(/2).. .7r(/Ar)|0) (6.45) 
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where the /’s are drawn from the same Gerfand-Shilov space Sa- 


The need to smear the correlation functions in order to give finite values to them in 
position space is an indication that the theory exhibits some degree of non-locality. The 
precise degree to which it is non-local depends on the precise GelTand-Shilov space needed 
to define these correlators. This in turn is determined by the precise exponential growth 
of the spectral densities. In other words, it is the exponential growth of the spectral 
density which is the signature of non-locality. By contrast, in a local field theory these test 
functions can be taken arbitrarily close to delta functions. 


6.5 Direct calculation of Galileon Spectral Density: (tt, x ) frame 

We now come to the computation of the spectral density. We shall do this for simplicity in 4 
dimensions (for which 04 = l/(27r^)), although the following argument is easily generalized 
to any dimensions. Our naive expression for the two point function 
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implies that the spectral density is given by 
e{k°)9{-k‘^)27rp^{-k‘^) = j d^ F{y,k) 
Where, by Ijy"^ we mean 


exp 


27r^A® \yt. 2/® 




[27r^A® 


y- 


(6.46) 


(6.47) 


y2n _ ^yO _ 

To make sense of the integral we perform an expansion in the coupling constant 1/A® 


(6.48) 


9{k^)0{—k‘^)2TTpT^{—k‘^) 


dSe-^^-yF{y,k)Y, 

n=0 


1 

n\ 


1 

27r^A® 


fk^ (k-y)^ 

W- y^- 


n 


(6.49) 

The validity of this expansion will be justified after the fact by demonstrating that the 
spectral density is an entire function of 1/A®. Neglecting the contribution from F{y,k), 
then each term in the sum is associated with contributions from n particle intermediate 
states. This is consistent with the definition of the Kallen-Lehmann spectral density, each 
term must be positive. Now F(y, k) contains a finite number of terms of which the first 
one is l/(47r^)yZ^. Let us consider this first contribution 


r 1 

9{k'^)9{-k‘^)2TTp^fl{-k^) = J ^ 


1 

n\ 


1 (k-yf 

A ^ 


= / d%e-^^-y 


4:Tr‘^y_ 


00 ^ 

T- 

„ n! 


2 

n=0 
00 n 


1 


27r2A® 


n=0 

^_k^) ^ ^^.^k ■ yf 


2'k‘^K^ \yt. 


y- 


y- 


y- 




(-!)’■ 4*^ 


{-k^Y{-ik.y) 


47 r 2 fUn — r)!(27r2A®)’' 

n=0 r=0 


2n-2r 


^6n-2r+2 


n 


(6.50) 

(6.51) 


- 57 - 



Now by a simple extension of the result in Appendix A 
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we may by repeated differentiation with respect to a followed by setting a = 1 infer that 
'' {-ik.yf^ 


and so the spectral density is given by the double sum (using N = 3n — r + 1), 
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which at large /r is approximated as 
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The full contribution to p coming from the additional terms in F can be evaluated using 
the same method with a qualitatively similar result. In particular the full p{p) is manifestly 
positive guaranteeing unitarity, and has order of growth a = 3/5. This same method may 
then be easily generalized to any dimension to give a = {d + 2)/{2{d + 1 )). 


As we have seen, the spectral density is an entire function of exponential growth. One 
way to understand the exponential growth is to use the fact that at large p the sum will 
be dominated by large n, and so may be safely approximated by an integral. The growth 
of the term is compensated by the fall off of the factorial coefficients in such a way 
that there is a peak value of n which determines the leading contribution to the sum. 
Performing a saddle-point approximation the peak value of n is found to be 


n ~ 



(6.57) 


However, the sum over n corresponds to the sum over n = A^-particle intermediate states 
in the Kallen-Lehmann spectral density. We thus see that the resonant contribution to the 
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spectral density of the Galileon operator comes principally from A^-particle intermediate 
states (generalizing to d dimensions) with 

f E \ 

-/V ~ ^ j ^Er^{E). (6.58) 

with E ~ The typical energy of a quanta is E/N = l/r^{E), i.e the frequency associ¬ 
ated with the Vainshtein/classicalization radius. 


Although this is a free theory, and there will be no scattering, assuming this property 
to hold true for the Wightman functions of an interacting theory, and given the relation 
between the scattering amplitudes and these functions, it is reasonable to suspect that the 
scattering amplitudes for scattering of two hard quanta with energy E will be similarly 
dominated by iV-particle out states with N ~ Er^,{E). This is precisely the conjecture of 
the classicalization scenario [23-26] [46] . There the intermediate states are viewed as bound 
states of N soft quanta which ultimately decay into free quanta. These bound states may 
be viewed at large energies as approximately coherent states, i.e. solutions of the semi- 
classical equations of motion. This is precisely what we have above, except the states are 
not truly bound since we are dealing with a free theory. 


To reiterate the central point, we have confirmed the exponential growth of the spectral 
density, which implies that there is a fundamental bound on the locality of the theory: local 
effective field theory is recovered only for 

1 /E\ (sii) 

|x|>r*(.E)=A f—j . (6.59) 

Nevertheless, unitarity and Lorentz invariance are maintained up to arbitrarily high ener¬ 
gies. 


6.6 Direct calculation of Galileon Spectral Density: {p, x) frame 

In the dual (p, y) frame, the dual Galileon field p is simply a free field with a standard 
two-point function. On the other hand, the matter field y satisfies a non-trivial equation 
which is dependent on p. The non-localizability does not disappear in this frame, but is 
contained entirely in the matter field spectral densities. To compute this we use the formal 
operator definition given in Sec. 6.3 from which we infer 
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Then defining 


we have 
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Performing the integral over k' we have 
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The spectral density is then 


= 2.S(,)+± i (^)" • 

(6.70) 

In particular in four dimensions we have 

27rp^{p) = 27rS{p) + ^ • (6-71) 

This result is derived from a different point of view in Appendix C. Explicitly this is 
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once again confirming the predicted growth a = 3/5. This demonstrates that regardless 
of which duality frame we use, the field theory will remain non-localizable with the same 
order of growth when there are more than one interacting species. 
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6.7 Prom Wightman to Feynman to Schwinger 

The calculation of the S-matrix requires knowledge of the time ordered correlation func¬ 
tions (Feynman) or equivalently the generalized retarded products. Once the Wightman 
functions are known then the time ordered correlators determine the retarded products 
uniquely (and vice versa). However both require an off-shell extension not needed in 
specifying the Wightman functions. The existence of well-dehned Wightman functions in 
momentum space is not sufficient to uniquely specify the Feynman or retarded correlators. 
That is because their usual dehnition requires a Lorentz invariant notion of time-ordering. 
However, in the case of a non-localizable field theory the position space correlators do not 
exist and so there is no sense in which we can time order them to define the time-ordered 
or retarded products. 

There is another way to see this problem, which also gives us a clue to the solution. 
Consider an ordinary local field theory for which the Wightman function spectral density 
growth is bounded polynomially by p < for integer N. A naive definition of the 

Feynman propagator in momentum space would be given by 

poo _■ 

This is a dispersion relation which shows that the Feynman propagator is an analytic 
function of z = — in the whole complex z plane, modulo a branch cut along the real 
z > 0 axis (see Fig. 2). 

poo ■ 

Gf{z)= / d|U/9(^) -. (6.75) 

Jq z — /i 

with the physical propagator defined as the limit to the real axis from the upper half z plane. 


However this expression is ill-defined as it diverges at large ^ due to the polynomial 
growth of the spectral density. This divergence shows up in the imaginary part only, the 
real part, being determined by the pole at /U = —is finite. As is usual for dispersion 
relations the resolution is to perform a finite number of subtractions at some reference 
point = piQ 
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(6.76) 


where the coefficients Cn{po) are formally defined by 

= (6,77) 

The divergence of the original expression is now confined to these finite number of coef¬ 
ficients and they may consequently be renormalized by subtracting infinite counterterms 
which amounts to allowing the Cn to be arbitrary finite coefficients. 
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Figure 2. Contour used to define Dispersion relation for Feynman propagator. The branch cut 
and pole are shown in red. 


In position space these terms correspond to derivatives of delta functions, so-called 
contact terms, as they are all localized at the coincidence or ‘contact’ point x = y. They 
are thus all associated with local counterterms needed to subtract UV divergences. The 
implication is that even in a local field theory whose spectral density grows only polynomi- 
ally, the time ordered correlation function for two fields is defined only up to counterterms 
which renormalize the contact terms even once the Wightman function has been deter¬ 
mined. 

The freedom of choice of these renormalized contact terms is the same ambiguity that 
resides in defining a time-ordered product. For instance, the time-ordered product between 
two fields A and B can be expressed as 

a2n 

TA{t)B{t') = e{t - t')A{t)B{A) + - t)B{t')A{t) + Y, Cn{t)^5{t - t'). (6.78) 

n 

Since all the contact terms, vanish away from t — t' = 0, any expression of this form is 
an acceptable definition of a time-ordered product. An explicit example of the ambiguity 
in defining time ordered products is the difference between the normal T product used 
in canonical Hamiltonian perturbation theory and the T* product used in the covariant 
Lagrangian perturbation theory. The equivalence of the Hamiltonian and Lagrangian def¬ 
initions of perturbation theory rests precisely on the existence of those nonzero contact 
terms which arise in the different definitions of time ordering (this is sometimes known 
as Matthews theorem [124], a more general discussion can be found in the textbook of 
Bogoliubov and Shirkov). 
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An equivalent way to express this ambiguity is to introduce a polynomial indica- 
tor/indicatrix function g{—k‘^) = which has the same polynomial growth 

at high k as the spectral density and has no zeros for positive real arguments dn > 0. Then 
we may define the Feynman propagator via 

Gpik) = g{-k^) [ dg — -+iu{-k^). (6.79) 

Jo 9\k-) k^ g — le 

where u{—k‘^) is a real polynomial in kJ. To see the equivalence of this with the previous 
expression it is sufficient to note that 
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k^ + g — ie 




n=0 


(fc" - goT 
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g{g){k‘^ + g - ie) ’ 


(6.80) 


where g{g) = {g + go)^ ■ The polynomial u{—k‘^) accounts for the contact terms. More 
generally it is sufficient to take g{g) to be any A^-th order polynomial whose zeros do not 
lie on the positive real axis. Reversing, the argument we see that for a generic polynomial 
Ai-th order polynomial g{g) we can write 


k^ + g-ie g{g){k‘^ + g - ie) ^(n + 1)! g{g) 


(6.81) 


Modifying the form of g{g) just modifies the contact terms and so any real polynomial 
g{g) is an acceptable choice. 


It is straightforward to see that regardless of the precise form of g{g), as long as g{g) 
is real the on-shell part is 

Re[GF{k)] = Tip{-k^). (6.82) 

The off-shell part affects only the contact terms. We may see this easily by passing to 
position space for which the Feynman propagator will take the form 

Gf{x) = g{D) {e{x^)Wg{x) + e{-x^)W*{x)) + iu{n)5*{x) (6.83) 

where we have dehned the regulated Wightman function 

tU(U = / ■ (6-84) 

Due to the indicator function, Wg is well-defined in position space. Since it is Lorentz 
invariant we may time order it in the standard way. Then the polynomial differential 
operator g'(n) acts on this time ordered Wightman function giving two sets of terms. The 
action of g{a) on Wg will just give back the original Wightman function, whereas the 
action of g{0) on the 6{x^) functions will give contact terms. Since <?(□) is a polynomial 
there are a finite number of such terms and so we find 

Gf{x) = [6{x^)W{x) + 0{—x^)W*(x)) -|-hnite number of contact terms. (6.85) 
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Generalization to non-localizable field theories 

In passing from a standard local field theory to a strictly localizable or non-localizable 
field theory, we pass from needing a finite number of contact terms (counterterms) to an 
infinite number. Mathematically this corresponds to replacing the polynomial g{—k‘^) with 
an entire function that satisfies the same requirements that has the same order of growth 
as /o(/i) and has no zeros for —real and positive. Then we may define the universal form 
of the momentum space Feynman propagator as 

Gpik) = gi-k"^) [ d/i 1,2 ^ — -+iu{-k‘^), (6.86) 

Jo 9W k^ + n-te 

where u{—k‘^) is now an entire function. Once more this definition is unique on-shell, and 
all the ambiguity is contained in the off-shell part. It is naturally extended to the complex 
plane z = —k‘^ via the dispersion relation 

Gf{z) = g{z) [ — - - \-iu{z). (6.87) 

Jo 9W z- 

In the case of non-localizable field theories this expression cannot be directly Fourier trans¬ 
formed without appropriate smearing with test-functions. However, in momentum space 
it is well-defined. 


The analogous definition of the retarded propagator is now uniquely specified in terms 
of the Wightman function and Feynman propagator to be 


Gnik) 


GF{k)-W{k)=g{-k^) 



pjk) _ zi _ 

9{p) kJ — {kP + ieY + g. 


-|- iu{—k^). 


( 6 . 88 ) 


Since the function g{—k‘^) is entire both the Feynman and retarded propagators have the 
same analyticity properties as for a standard local field theory. However, since they are 
not polynomial bounded we may not give the standard finite number of subtractions dis¬ 
persion relations for them. Since it is apparent that different choices of g{—k‘^) will give 
rise to different u{—k‘^), we may set u{—k‘^) to zero, understanding that this extra freedom 
is already built into g{—k‘^). 


This procedure may be generalized to give time-ordered correlation functions for any 
number of fields [125, 126]. Consider a N-point Wightman function which may be expressed 
in momentum space in the form 


(0|(^(xi) . . .(^(xAr)|0) 



Wiki ...%). 

.n=l / 


(6.89) 


We introduce a symmetric N-point indicator function g{ki,... k^) which is an entire func¬ 
tion in all its arguments, has no zeros in the physical region, and has the same growth at 
large arguments as W{ki ... /cat), then the time-ordered correlation functions in position 
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space are given formally by 


Gf{xi, .. .,xn) = - idN) 


and 


0{xi - X 2 ) ... 0{XN-1 - xn) Wg{xi,..., xn) 

permutations 

(6.90) 


Wg{xi, .. .,xn) 


d^xi ... 


d^XAT kn.Xn^ii) 



W{h...,kN) 

g{ki, ...,kN) 


(6.91) 


To these we add entire function contributions which vanish on-shell. As usual these formal 
expressions become well-defined when appropriately smeared with test-functions. Here the 
sum over permutations is the standard time-ordering operation. The anti-time ordered 
correlation functions may be defined analogously and are given by the complex conjugates 
G*p{xi,... ,xn). The generalized retarded products which are usually formally defined as 

'jz{x, XI, . . . , X„) = (-i)” ^ ^(x-Xl) . . . 6'(x„_l-Xn)[. . . [[(i>{x), ${xi)], ^{X2)] . . . ${Xn)] , 

Permutations 

(6.92) 

are determined unambiguously once the time ordered products have been specified 


7t(x,xi,... ,Xn) = 
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G*p{xk, ..., Xn)GF{x, Xi, . . . , Xk) 


(6.93) 


Schwinger functions 

Having constructed the time-ordered correlation functions, and using the fact that the in¬ 
dicator functions are entire and hence well-defined under Wick rotation, it is now straight¬ 
forward to Wick rotate to give the Euclidean correlation functions (Schwinger functions) 


Ge{xi, ... xn) = - idN) [GE,g{xi,... xat)] , (6.94) 

where GE,g{xi,.. .xn) is the Wick rotated version of 


0{X1 - X2) . . . ^(xAr-i - xn) Wg{xi, . . . xn) 

permutations 


^Wick rotation^ GiJ,3(xi, . . . Xjv) • 

(6.95) 


In particular, in the case of the two-point function, we have in momentum space 

p{k) 1 




(6.96) 


This is the obvious Wick rotation of the Feynman propagator in momentum space. It is 
not, however, the analytic continuation of the Wightman function since the Wightman 
function contains no reference to g{—k‘^). This is an important difference with standard 
local field theories, the latter having the property that the Wightman function at space¬ 
like separations is the Euclidean correlation function. We also see that it is impossible 
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to uniquely specify the Schwinger functions without reference to the indicator. This is 
a reflection of the fact that Euclidean theories are always off-shell, and are consequently 
always sensitive to the ambiguity of the off-shell extension. In path integral language this 
is related to the inherent ambiguity in specifying the measure of the path integral. 

The most striking consequence of this distinction between the Schwinger function and 
the spacelike Wightman function is that for the case of a non-localizable field theory, 
even though the spacelike Wightman function may not exist in position space without 
appropriate smoothing, the Euclidean correlator can exist since it is determined by the 
behaviour of g{—k‘^) for k'^ > 0, i.e. for negative arguments. Even though g{—k‘^) grows 
exponentially for large timelike k, for certain theories it may decrease exponentially or as 
a power for large spacelike k. An example of a non-localizable spectral density with this 
property is p{g) = for which we may choose g{—k‘^) = . This renders the UV 

behavior of the Euclidean theory much better than a standard local field theory. This 
fact is the basis of Efimov’s non-local field theory program (see e.g. [86, 127]) which 
utilizes this improved UV behavior to regulate loops in a perturbative expansion. Our 
perspective is somewhat different though since we imagine that this behavior comes from 
the nonperturbative properties of a classically local field theory and it is not clear that 
there is any sense in which we may find a perturbative expansion in which the Euclidean 
propagator is replaced by an expression of the above form. More generally, there is no 
expectation that g{—k‘^) decreases exponentially at spacelike /c^. It may, for example, 
simply oscillate. 

7 Superluminality, Macro-causality 

7.1 Causality and the Jost-Lehmann-Dyson representation 

One of the most commonly used arguments against the possibility of UV completing 
Galileon theories and massive gravity theories in a Lorentz invariant manner is the ap¬ 
parent existence of superluminal solutions in the low energy effective field theory. Indeed 
the superluminal solutions that arise in massive gravity are already present in its Galileon 
decoupling limit (for a discussion on this connection see the review [7]), and it is sufficient 
to analyze the decoupling limit theory to understand these. 

These arguments can be stated in slightly different but entirely equivalent ways 

• The LEEFT admits classical solutions whose perturbations have superluminal group 
and/or phase velocities. 

• A characteristic analysis of the classical equations of motion implies that there are 
initial data for which the Gauchy problem becomes ill-defined - the ‘Velo-Zwanziger 
problem’. [128] (For a discussion of this for massive gravity see [129, 130]). 

• There exist solutions for which the momentum-velocity relation for the field may not 
be inverted. As such there is a conflict between the desire to canonically quantize the 
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theory using equal time commutation relations and Lorentz invariance - the ‘ Johnson- 
Sudarshan problem’. [131] 

As a consequence, these results are often argued to be fatal to any attempt to UV complete 
such theories, at least in a Lorentz invariant way. It is argued that: 

• Lorentz invariance is fundamentally broken in these theories. 

• The existence of superluminal propagation implies the formation of Closed Time-like 
curves and hence the violation of causality. 

Although all of these statements are correct classically, their importance in the quan¬ 
tum theory is more subtle. To begin with, neither the low energy phase nor group velocity 
has anything to do with causality, rather as is well known this is determined by the front 
velocity which is the high momenta limit of the phase velocity vp = limfc_,.o up(A:). The 
characteristic analysis correctly determines the classical front velocity, however since this 
velocity is determined by the high frequency limit, in theories built on the Vainshtein 
mechanism the field is highly strongly coupled there, and the tree level calculation cannot 
be trusted (we shall see this explicitly below). The fact that the momentum-velocity re¬ 
lation cannot be inverted is not a problem semi-classically since those solutions for which 
it cannot be inverted will not be described by the LEEFT. Quantum mechanically, this 
observation is consistent with our earlier arguments (Sec. 4.4) that Vainshtein type theories 
should not be canonically quantized and, as we have emphasized throughout, there is no 
obstruction to quantizing in a Lorentz invariant manner. Finally, the classical formation 
of closed time-like curves is never in the regime of validity of the LEEFT theory [132]. 


A concrete realization of the apparent superluminal propagation can already be seen 
in the case of the special Galileon which is dual to a free theory. In any dimension, the 
dual Galileon theory exhibits classical solutions which exhibit superluminal propagation 
(for details on these arguments see [74, 76, 77, 118]). In our current language this is the 
statement that if we evaluate the two point function of the operator 7r(x) in a coherent 
state [a), and maintain only the contributions that arise at tree-level, then there will exist 
some coherent state |a) for which the tree-level two point function has a commutator which 
is non-zero outside of the Lorentz invariant light cone. In fact, the simplest example of this 
is in the (p, x) dual frame. Evaluating the two-point function for x in a coherent state for 
p we obtain 


(a|x(a;)x(2/)|a) 


J dk exp 
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(7.1) 


where 

{a\p{x)\a) = pc{x) 


(7.2) 


defines the background field in the (p, y) dual frame. This clearly propagates in a lightcone 
determined by 
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(7.3) 
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which can easily be chosen to lie outside the usual lightcone (x — y)'^ < 0 for some choice 
of Pc{x). This example is discussed in more detail in [118]. However, the above is only a 
tree-level calculation, and the all-loop (orders in h) order answer 


{a\x{x)x{y)\a) = 


j dk exp 


ik. 



dpcjx) 


dPc{y) \ 
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exp 



(fc- {x-y)f \ 

(a;_y)d+2 j 


(7.4) 


is only meaningful after smearing with test functions. The additional exponential which 
arise from loop corrections (we have included a factor of h to make clear its dependence) 
completely modifies the high k behaviour of the correlator, and this is relevant since causal¬ 
ity is determined by the front velocity which is the high k limit. The above expression shows 
clearly that the tree-level approximation cannot be trusted in calculating this. 


To analyze this two-point function, as in the vacuum case, it is helpful to work with 
a spectral density which is finite. Since we are no longer in vacuum, the relevant spectral 
density is the one provided by the Jost-Lehmann-Dyson (JLD) representation [133, 134]. 
The relevant physical question is: Does the UV completion of Galileon theories 
allow for macroscopic superluminal propagation? Here, by macroscopic, we mean 
at distances larger than the locality bound. This distinction is important since in a trivial 
sense superluminality occurs in the region forbidden by the locality bound since their is no 
notion of locality there, but this superluminality is fundamentally harmless. Provided that 
the theory respects a Lorentz invariant notion of macroscopic causality, then there will be 
no real superluminal propagation. The tree-level calculation on the other hand implies a 
superluminal propagation which would be seen at arbitrary large distances. 


A meaningful definition of macroscopic causality is as follows. We define a pair of test 
functions fxoix), gxoix) from the space Sa which are localized at the point x = xq to the 
maximum extent allowed by the locality bound of the theory. We define the smeared field 

Hfxo) = j (i'^xfxo{x)$ix) . (7.5) 


Then macro-causality requires that 


Hfxo): ^{9yo) 0 > (^0 - yof +00 


(7.6) 


for two such localized test functions /, g, i.e. the vanishing of the smeared commutator 
at large space-like separations. Furthermore this vanishing should be exponentially rapid. 
We further require this for all allowed pairs of test functions from the space Sa, including 
those test functions obtained by a Lorentz boost from the original one. In this sense, this 
definition of macro-causality is Lorentz invariant. This condition alone is sufficient to for¬ 
bid superluminal propagation since any superluminal propagation would inevitable lead to 
an enlargening of the causal cone at large distances, leading to a contribution which falls 
off only as a power law. 
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Although superluminality is usually diagnosed around classical backgrounds, i.e. in 
coherent states, it is simpler at the quantum level to work with momentum eigenstates. A 
momentum eigenstate with a finite number of particles will not exhibit any superluminality, 
and so we must still work with states with infinite number of particles as in the case 
of coherent states. This is easily achieved by projecting the coherent state \a) onto its 
momentum eigenstate as follows 

\a,P) = I d^ye-^^-y\a'y), (7.7) 


where a'y is the translated coherent state whose vev is (j)c{x + y) if a has vev 4>cix), so that 
|q;,P) defined a momentum eigenstate containing an infinite number of particles. More 
precisely, if the original coherent state takes the normalized form 


then 


1^) = g-UdfcKfc)Pexp 

\a,P)= [ 


dk a{k)al. 


| 0 ), 


exp 


dk a{k)P^’yd\ 


| 0 ). 


The original coherent state is obtained as a superposition of such states 

f d'^P 


a = 


\a,P) . 


(7.8) 

(7.9) 

(7.10) 


. (27r)'^' 

Given this, to address causality it is sufficient to analyze the commutator between two such 
coherent momentum eigenstates. Denoting these states by the further shorthand \Pi) and 
|P/) then we consider 

(P/l [m.m] \P,) = [^{z/2)A{-z/2)\ IP,), (7.11) 

where z = {x — y). 


In a strictly localizable field theory we would require that (P/| |^())( 2 :/ 2 ), (/>(— 2 ;/ 2 )J |Pj) 
vanishes outside the lightcone < 0 for all |Pj) and |P/). This would be sufficient to ensure 
the absence of superluminal propagation and micro-causality. In a non-localizable theory 
we demand the weaker condition of macro-causality. A sufficient condition for Lorentz 
invariant macro-causality is that 


(P/I ^izl2)J{-zl2) \P,)=g{a,)F{z,Pi,Pf), 


(7.12) 


where F{z, Pi, Pf) vanishes outside the lightcone < 0 and g{fi) = h{g)^ is an entire 
function of the desired order of growth built out of the square of an entire function h{g). 

To see why this is sufficient consider the smeared commutator 

(P/I [^{UA{gy,)] IP) = [^{zl2)A{-z/2)\ |P) 

= y d'^'xy d!^zU,{x + z/2)gy,{z - z/2)P^P^-Pf)-h{U,)h{U,)F{z,Pi,Pf), 

= I d'^x I d^zfl,{x + z/2)gy,{z-z/2)P^^^-^f^^F{z,Pi,Pf), (7.13) 
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with X = {x + y)/2 and we have defined the new test functions 


(7.14) 

(7.15) 


/^(,(x + z/2) = h{a^)f^g{x + zl2 ), 

Qxoix - z/2) = h{n^)g^^{x - z/2). 

If the original test functions are drawn from the space 5 q,, then the new test functions 
will remain localized at xq and yo respectively (i.e. the growth of the Fourier transform of 
h{k‘^) is compensated by the fall of of the Fourier transform of the test function). 

Given this as (xq — yo)^ = —>■ +oo, the product of test functions fxQ{x + z/2)gyQ{z — 
z/2) will drop off exponentially away from z"^ = L?' ^ +oo. However, by our assumptions 
F{z,Pi,Pf) = 0 for > 0. Consequently, the only contribution to the integral 

comes from the tail of the overlap of the test functions which tends to zero exponentially 
as (xo — yo)^ = —>■ +oo, thus fulfilling the requirements of macro-causality. 


Thus to demonstrate Lorentz invariant macro-causality for a given field 0(x), it is 
sufficient to demonstrate that the following quantity vanishes outside the lightcone: 


F{z,Pi,Pf) 


2 1 
' W) 


d^y - 


( 2 vr) 


ik.y 


{Pf\ <^(y/2),<^(-y/2) 


|Pi)=0, z^<0. 


(7.16) 

To prove that this is the case, let us consider the (tt, x ) frame, and consider the local 
field 7r7v(x) defined in Eq. (6.15). Since this field is local, and its correlation functions are 
tempered distributions, then its commutator vanishes outside the lightcone 


[7rAr(x),7rAr(y)] = 0, (x-y)^<0. 


(7.17) 


Given this, it is subject to all the requirements for there to exist a JLD representation 
[133, 134]. We may thus write 

roo r- jdu f 

{Pf\[nN{y/2),TrNi-y/2)]\P/) = J dy J j P^uP^-yDN{u,y,Pi,Pf)2TTe{k^-u^)5{{k-u)'^+g) 

(7.18) 

where y, Pi, Pf) is the JLD spectral representation which depends on the specific 

momentum eigenstates and is the analogue of the Kallen-Lehmann spectral density when 
not in vacuum. Since tttv is a tempered field Dn{u, y, Pi, Pf) grows at most as a polynomial 
in y. This is sufficient to ensure the convergence of the integral since 

(P/|[TO{!,/2),iA,(-!,/2)]|F.) = j ^e“-»0«(<.,-(i-<.)ta,P/)2>rt(t“- «"), 


for which the k integral is convergent. Hence for this field we have 

P^k f d<^u ^,,,DN{u,y,Pi,Pf)g{-{k-u)^) 


FN{z,Pi,Pf) = / dy 
Jo 


{2tt)<^ j ( 27 r )“ 


gi-k"^) 


gik-) 


2Tre{k'^—u^)5{{k—u)‘^+y) 

(7.20) 
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Now this can be rearranged in the form 


FN{z,Pi,Pf) = 


dfj, 


iu.zF>N{u,^,Pi,Pf) g{-{-id^f) 


10 


(2vr) 




g{-i-idz + u^) 


A^{z), (7.21) 


where A^{z) = f -^^e^^'^2'ne{k^)5{k‘^ + g) is the free-held of mass commutator which 
automatically vanishes for < 0. 


We may now take the limit N ^ oo, since although 


D{u, g, Pi, Pf) = lim Dn{u, g, Pi, Pj), 

N^OO 


(7.22) 


will have an exponential growth in this limit D{u, g. Pi, Pf) ~ this will be compensated 
by the indicator function g{g). As we have already discussed, any expression which remains 
hnite in the limit N ^ oo will give the appropriate dehnition for the Galileon held 7r(x). 
Thus we require that the following expression is hnite and vanishes outside the lightcone 


roo p 

F{z,P„Pf) = j^ dgj 
First let us consider the part 


iu.z F{u,g,Pi,Pf) g{-{-id^f) 
{2 -kY^ g{g) g{-{-idz + u)"^) 




9(-Hd.r) 

g{-{-idz + uy) 



d'^k g{—k‘^) 
(27r)^ g{—{k + u)^) 


27re(fe°)<5(fe2 + . 


(7.23) 


(7.24) 


Since g{—k‘^) is an entire function, then the ratio is an analytic function over 

the whole complex plane modulo poles at gi. We may thus write 


g{-{k + u)^) ^ -{k + u)^-gi' 


(7.25) 


where h{k, u) is an entire function of k and Ri{k, u) are the residue terms at the poles which 
are also entire functions of k. Here we have assumed that all of the poles are simple poles 
for simplicitly (we can always view multiple poles as limits of combinations of simple poles). 


In any contour integral, the contributions from the poles are harmless since their 
residues will give rise to contributions to F{z, Pi, Pf) that vanish when acted on by g{—n,^). 
This is just a rehection of the fact that the only reason these poles exist is because we have 
chosen to extract the operator g(—Dz) outside of the integral, but the inverse of this 
operator is ambiguous due to the poles. 

To make this more precise we note that we can write 


g{a)F{z, Pi, Pf) = g{D)F,{z, R, Pf) + g'i{a)F^{z, R, Pf) 


(7.26) 


where g'i{g) = g{g)/{g — gi) is an entire function with a given zero removed, and 

d"'^ iu.zF{u,g,Pi,Pf) 


Fi{z,Pi,Pf)= / d^ 
Jo 


Fi{z,Pi,Pf) = 


dg 


{27r)^ g{g) 

d'^^^ iu.zF{u,g,Pi,Pf) 


(2vr) 


9{9) 


h{-idz,u) A^{z), (7.27) 

Ri{-idz,u) A^{z). (7.28) 
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In the case of higher order poles of multiplicity N a similar approach may be taken in 
which = g{g )/{^l — and Ri{k,u) is replaced by the appropriate higher order pole 
residue. We may now refine our statement of macro-causality that independently we must 
have Fi{z, Pi, Pf) and F^iz, Pi, Pf) vanish for z^ > 0. 


The key observation is that, if g{fj,) is an entire function of // of order a, i.e. g{g) ~ 
then both h{k,g) and Ri{k,g) are entire functions of k of order 2a — 1. To see this, it is 
sufficient to look at the fixed u, |/c| —)• oo asymptotics of 


g{-{kPuY) 


^ g<7(-fc2)-_^(_(fc+n)2)“ ^ ^2aa{-Pr-^k.u < ^2aa\k\^-‘-^\u\ 


(7.29) 


Thus, if g{—k‘^) has the exponential growth of a non-localizable field with 1/2 < a < 1 then 
the ratio has the exponential growth of a strictly localizable field. For instance, 

in the case of Galileons in four dimensions, a = 3/5 and so 


~ , (7.30) 

g[-{k + uY) 


and hence both h{k, u) and Ri{k, u) are independently bounded by However, we 

know any entire function of —id acting on a function which vanishes outside the lightcone 
and has a growth less than a linear exponential in k, will still vanish outside the lightcone 
because that is how we defined strictly localizable helds in the first place! [78-80]. 


To see this more precisely, consider the action of h{—id, u) on the retarded propagator 

= 9in.z)A^{z ), (7.31) 

where is a future pointing timelike vector. As long as we can show that h{—id, tt)Gret,^('Z) 
vanishes for n.z < 0, for any future time-like vector n, then this proves the vanishing of 
h{—id, u)Gret,iJ,iz) outside the whole lightcone since this is the only way we have a Lorentz 
invariant (i.e. independent of n) notion of time ordering. The most dangerous term then 
comes from when all the time derivatives act on the 9{n.z) function since this gives the 
contribution most likely to extend outside the lightcone. Thus 

h{-id, u)Gret,tJ.{z) = {h{-id, u)9{n.z))Afj_{z) + ..., (7.32) 

However, viewed as a function of the time-like direction n.k alone, we know that h{k, u) 
is an entire function which grows at most as which being strictly localizable, 

is insufficient to give support for n.z < 0 [78-80]. Hence, since h{—id,u)Gret,fi{z) = 0 for 
n.z < 0 for all future time-like n, then h{—id,u)Gret,iJ.{z) = 0 for > 0 and by extension 
h{—id,u)A^{z) = h{—id,u) {Gret,iJ,{z) — Gret,iJ,{—z)) = 0 for > 0. Identically the same 
argument applies to Ri{—id,u)A^{z) = 0. 

To hnally show that Fi{z, Pi, Pf) and F^iz, Pi, Pf) vanish for > 0, it is sufficient to 
ensure that the remaining integrals converge. The g integral clearly converges since our 
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assumption is that the indicator function g(fj,) is chosen to have a growth to compensate 
the exponential growth of D(u, n, Pi, Pf). Thus, the only remaining concern is for the 
convergence of the u integral. Previously we considered the action of h{—id,u)A^{z) at 
fixed u. However, it is possible that h{k, u) could grow exponentially in u as h{k, u) ~ 
at fixed k. At the very least we have the following bounds 

h{—id,u)A^{z) < , (7.33) 

u)A^{z) < Bi{z)e^\^\"^ . (7.34) 

Thus a sufficient (but perhaps not necessary) condition for the convergence of the remaining 
integrals is that 

D{u,ix,Pi,Pf)e^\^\"\ (7.35) 

is polynomially bounded in u, or at least is bounded by a linear exponential growth in u: 

D{u, fi, Pi, . (7.36) 


In fact we may now withdraw from our use of momentum eigenstates and consider the 
commutator between coherent states |a) and |/3) given formally as 


(/I| [7r(x), 7r(y)]|a) = J J J e*”'^£)(u,/r, a,/I, x) 27re(A:^ — u^)(5((A: — u)"^ + ^u), 

° (7.37) 

where 


D{u, n, a, f3, x) = 


d^Pf f d^'P, 


(27r)'^ J {2ttY 


y^Pi-Pf)-D{u, g, (/3, Pf), {a. Pi)), (7.38) 


and D{u, iJ., (/?, Pf), {a. Pi)) is the JLD spectral function associated with \a, Pi) and |/3, Pf). 
Hence for a given initial and final coherent state, \a) and |/3), to exhibit no superluminali- 
ties, a sufficient (although possibly not necessary) condition is that 

b{u, fi, a, p, x)e'"l“l"“ < , (7.39) 


is bounded by a linear exponential in u. 


This is a restriction on the state itself since u essentially encodes the momentum of 
the background coherent fields. To see this we note that we can write 

r°° - 7 d'^A- 

{l3\['k{x),Tr{y)]\a) = J dgD{iJ,,a, P,x,y) J P^-^2TT€{k^)6{k'^ + y), (7.40) 

where 

r d^v 

D{y,a,P,x,y) = J y,a, P,x). (7.41) 

From this representation it is clear that D{y,a, P,x,y) is a local version of the Kallen- 
Lehmann spectral density, whose x and y dependence will arise only because we are not 
evaluating this expression in vacuum, but rather between two coherent states. Because of 
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this we do not expect D{ii, a, (3, x, y) to vary more in space than the background coherent 
state itself. In other words if {a\n{x)\/3) varies over a typical length scale Lq, then we 
expect the same to be true for f3,x,y) which in turn implies that the support of 

D{u, //, a, /?, x) will be for momenta with u ~ 

Now in the present case, the spatial dependence of D{fi, a, j3, x, y) will be determined by 
{/3\n{x)\a) which we know to be a solution of the classical ((i+l)’th order Galileon equation 
of motion. Alternatively we may use the duality to infer (/3|7r(x)|a) from {f3\p{x)\a). 
Consider a coherent state of total energy E. Although the Fourier transform of {f3\p{x)\a) 
may fall off only polynomially, we expect that the Fourier transform of {l5\Tr{x)\a) to fall 
off exponentially for \k\ ^ l/r*(Fl). The reason is that the classical equation becomes 
strongly interacting at distances |x| < r.^{E) which will suppress the momentum support 
for \k\ ;§> l/r^{E). Indeed following the classicalization arguments, we expect a state of 
energy E to be build out of N soft quanta with typical energy/momenta \k\ ~ E/N with 
N = Er^:{E). This behaviour will simply follow from the classical equations of motion. 
This in turn implies that the support of D(u, p, a, /?, x) is localized on modes u for which 
|tt| ~ E/N ~ l/r*(£'). As long as E ^ A then this support will be for w <C A and so 
we may reasonably expect that D{u, p, a, /3, x) decays exponentially for u ^ A. From the 
duality map 

m{x)\a) = I dS I ^U{{my)\o^)) (7.42) 

we may reasonably estimate using a saddle-point approximation that this exponential de¬ 
cay will be at least as fast as which is sufficient for convergence of all the 

integrals. We have not proven that all states satisfy this property, it nevertheless seems 
reasonable that generic coherent states fall off in this manner. 

We regard this as yet another manifestation of the locality bound. Those classical back¬ 
ground states that are not built out of high energy quanta localized in a region less that the 
associated locality bound will give rise to macroscopic causality. We have not (yet) rigor¬ 
ously excluded the existence of states which violate boundedness of D{u, p,a, (3,x)e^^'^^ 
and hence potentially give superluminal propagation, however such states attempt to local¬ 
ize quanta more than allowed by the locality bound. To put this more pragmatically, any 
semi-classical solution of the Galileon LEEFT that is built out of modes only with k < A 
should satisfy the condition that D{u, p, a, /3, is polynomially bounded. This is 

true even for configurations of total energy E ^ A since the power will effectively be cutoff 
in the strong coupling region \ k\ ^ l/r*(E) for which l/r*(E) < A. We will leave to future 
work a more precise exposition of this condition and its necessity. For now, we simply note 
that it is sufficient, and that there are an infinite class of non-trivial background solutions 
for which no macroscopic superluminality occurs in the fully quantum theory, even though 
the tree-level calculation would imply that there was. Any superluminality that does occur 
would appear to be intimately connected with the violation of locality due to the locality 
bound. 
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Extension to Retarded products 


The previous argument showed that quite generally the commutator of two fields vanishes 
at space-like separations. However causality requires an even sharper statement that we 
can define a retarded product of operators that has support only in the future light-cone. 
Our macroscopic version of this statement is that 

{Pf\R4'{fxo)^i9yo)\Pi} 0 > (xo- yof +00 , 

and (xo — yo)^ with — y°) < 0, (7.43) 


where R denotes the retarded product whose definition requires subtractions. The naive 
unsubtracted definition of the Fourier transform of the retarded product in terms of the 
JLD spectral density is 

/ /* OO /* J tZ 

d^ye-^^-y(Ffmy/2)<P(-y/2)lR} = dy j ^ P/) 

and so the subtracted form is 


— {kF' — -|- ieY + y 

(7.44) 


d<^ye-^>^-y{Pf\R${y/2)$i-y/2)\R 

rcsD 

/ dy 


/o 




—t 


= 9{-k ) dy 


{ 27 ry ^ g{y) (fc _ ^)2 _ (^0 _ ^0 + ^ 

d'^u g{-{k-u)‘^) D{u,y,Pi,Pf) -i 


(7.45) 


{2ttY g{-k^) 


5'(^) 


{k — u)2 — {k^ — + ie)^ -|- y 


up to the addition of entire subtraction functions. 

The argument proceeds as before: in order to satisfy macro-causality it is sufficient to 
have 

{Pf\R${y/2)$i-y/2)\P^ = y(n,)F^(y, R, Pf) (7.46) 

where Ffi{y, Pi, Pf) has support only in the future lightcone. Up to possible subtraction 
terms we have 


FR{y,Pi,Pf)= (7.47) 

r d^fc g{-{k - n)^) D{u, y, Pj, Pf) _ -i _ 

(27r)‘^ 7 (27r)^ g{—k‘^) g{y) _ {t )2 _ j'/j.o _ .^o _|_ ^^^2 ^ 

Following the previous reasoning, this is a linear superposition of free-space retarded prop¬ 
agators which will automatically vanish outside the future light-cone provided that the 
integrals converge and that ^ has the order of growth of a strictly-localizable field, 

which again requires only a < 1. The convergence of the u integral is ensured since 
de made polynomially bounded, and so the only doubt is the convergence 
of the u integral which will be satisfied under the previously stated conditions. 



- 75 - 



7.2 Causality and S-matrix Analyticity for Non-localizable theories 

A related statement of causality is that the S-matrix should be an analytic function of 
energy, modulo branch cuts associated with thresholds and poles associated with physical 
particles. We shall give a more thorough discussion of the properties of the scattering 
amplitude in [14], for now we give a simple example using the language of non-relativistic 
scattering that demonstrates our point of view that it is analyticity alone, and not poly¬ 
nomial boundedness, that is sufficient to preserve the notion of causality in a non-local 
theory. More precisely, we argue that analyticity alone implies macro-causality and only 
together with polynomial boundedness does it imply micro-causality. The latter condition 
is safely violated in a non-localizable theory, but the former can be maintained. We will 
give a more thorough discussion of these issues elsewhere [14]. 


To keep things as simple as possible let us consider a quantum mechanical S-matrix 
between two specihc states of energy E, Sfi{E), which is assumed to be analytic in the up¬ 
per half-plane, modulo the required poles and branch cuts, but non-polynomially bounded. 
Assume the existence of an entire function g{E) = ^2^ CnE^ which has the same order of 
growth at high energies as Sfi{E). Then, given our analyticity assumption, we may define 
a subtracted dispersion relation for the S-matrix element of the form 


/ OO 

-OO 


dE' 1 pfijE') 

27r g{E') E' -E-ie' 


(7.48) 


where Ufi{E) is an entire function of E. Ufi{E) may grow no faster than g{E) x polynomial 
so that the Fourier transform of Ufi{E) = Ufi{E)/g{E) is a tempered distribution. 


In the case of a local field theory it would be sufficient to take g{E) as a polynomial, 
then associated with this scattering matrix we could define a retarded operator 
dE 


GM=I 


—iEt 


SME) 


' —OO 

/■“ dE 


2tt 


/ C 

-C 


dE 
—( 
27r 


—iEt 


9 {E) 


dE' 1 


Pe{E') 


Loo 27r g{E') E' - E - ie 


(7.49) 


-(7.50) 


J—OO j—OO 

The first term, being the Fourier transform of a polynomial, can be expressed in terms 
of derivatives of delta functions. The second term can be organized into derivatives of 9 
functions. We then have 


GR{t) = Ufi 


d 


dt 


where 


5{t) - ig 
~ dE' 




-iE't 


. d 

'Et 

1 


HMt), 

(7.51) 

PfiiE '). 

(7.52) 


Hfi{t) manifestly vanishes for t < 0. Thus, again in a local field theory, Gr is built out 
of a finite number of derivatives of delta functions, which vanish away from t = 0 and a 
contribution which vanishes for t < 0 


GR{t) = Vfi 


.d' 

'dt 


6{t) — i9{t) f 
J —c 


dE' 

LtT 


-iE't 


PdE') 


(7.53) 
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where Vfi includes the contributions from the second term of (7.51) where at least one of 
the time derivatives acts on the 6{t) inside 

Thus we see that in a local theory, analyticity and polynomial boundedness would 
guarantee = 0 for t < 0. This is the usual statement of micro-causality as seen 

through the analyticity and polynomial boundedness of the S-matrix. This will continue 
to be true for strictly localizable field theories for which all of the above Fourier transforms 
are well-defined since the usual exponential oscillations from are sufficient to ensure 

convergence. Consequently we anticipate that the micro-causality condition holds for all 
strictly localizable theories 0 < a < 1/2. This has been proven for hxed physical momen¬ 
tum transfer in relativistic theories in [15]. 

In the case of a non-localizable field theory the micro-causality condition Gij(t) = 0 
for t < 0 can no longer hold due to non-polynomial boundedness of the S-matrix which 
is captured by the fact that the growth of g{E) ~ for some complex E grows too 

fast to ensure convergence of the contour integrals when a > 1/2. However, as we have 
seen, this also means this expression for does not exist, (at least not as a tempered 

distribution). This is quite a different situation than a genuinely acausal theory for which 
GR{t) does exist but has support for t < 0. 

However, what does exist in the non-localizable case is the causal response from any 
incoming state which is sufficiently delocalized in time, i.e. it is drawn from the same 
GelTand-Shilov space Sa as the test functions needed to dehne the operators. The formal 
response from that initial state is given by 

/ OO 

dt'G'ij(t - t')</'in(t') • (^-Sd) 

-OO 

or in Fourier space 

ME) = u„(£)*„(£) + g(E) ^ ^ e'’- E - 
We may rewrite this as 

= udE)UE) + , (7.56) 

where <pin{E) = g{E)(pin{E) and Ufi{E) = Ufi{E)/g{E), so that 

/ OO 

-OO 

Here G'^{t) is given by 


- 77 - 





Unlike GR{t), exists and vanishes for t < 0. This is the S-matrix equivalent of the 

causality condition imposed in Eq. (7.12). Micro-causality is only violated in the sense that 
J{t) has support in the past of 4>init), but only to the extent consistent with the locality 
bound. 

In this sense analyticity alone in the upper half-plane is enough to ensure macro¬ 
causality. We regard this is as sufficient to save us from the naive perils of closed timelike 
curves. Any closed timelike curves formed will be localized in the region where the locality 
bound forbids us from giving a well defined notion for time. Nevertheless macroscopically 
a well defined notion of time and time ordering emerges. This is consistent with the per¬ 
spective of Steinmann that analyticity in momentum/energy space is a sufficient notion 
for causality in a non-localizable theory [19]. This discussion should be contrasted with 
the case where analyticity is violated. For instance, imagine Sfi{E) has a simple pole in 
the upper half-plane at E = Eji + iEj. Then Gr will pick up an exponentially decaying 
contribution from the residue |G/j| ~ However this contribution blows up exponen¬ 

tially in the past destroying any notion of causality. A similar result would follow from the 
existence of branch cut (a branch cut may be viewed as a continuum of an infinite number 
of poles). 

7.3 An example of a non-localizable theory: Asymptotic Fragility of Fluctn- 
ating Infinitely Long Strings 

A simple example of a held theory which is consistent with our above reasoning is provided 
by the ‘asymptotic fragility’ for fluctuations in the effective held theory of inhnitely long 
strings [20] and, in particular, the model S-matrix considered in [21] which in our language 
is an explicit example of a S-matrix for a non-localizable held theory with a = 1. 

Consider the gauge hxed Nambu-Goto action for a string propagating in D dimensions 
in terms of its 77 — 2 transverse coordinates 

5 J d^a^-Det[rj,b - lldaX%Xi] . (7.59) 

In the case 77 = 24 this is of course the critical bosonic string. From the perspective of a 
two-dimensional held theory, this is, like the Galileon models, a non-renormalizable held 
theory for which in perturbative calculations we expect an inhnite number of counterterms 
to be included all at the scale Ig- Thus from the standard EFT logic we may only expect 
this theory to be meaningful for E <C 

Despite this, it is argued in [20] that the S-matrix for scattering of worldsheet pertur¬ 
bations on an inhnitely long string takes the remarkably simple exact form 

A{s) = , Im(s) > 0 , (7.60) 

A{s) = , Im(s) < 0 , (7.61) 

where s is the usual Mandelstam variable. This S-matrix, which is dehned by a single com¬ 
plex phase because we are scattering in two dimensions, is manifestly Lorentz invariant, 
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hermitian analytic ^*(s) = and crossing symmetric ^(—s) = A(s). It thus defines 

a UV completion valid at arbitrarily high energies E ^ and it can be shown at low 
energies to match the expectations based on the LEEFT (7.66). This example already has 
many similarities with our proposal for Galileons. In this case, the strong coupling scale 
A is 1/ls and the theory contains no UV fixed point, rather its high energy behavior is 
determined by l/lg- The authors of [20] refer to this as ‘asymptotic fragility’ to distinguish 
this behaviour from the UV hxed point of ‘asymptotic safety’. 


This high energy behaviour is precisely what is expected based on semi-classical scaling 
arguments. In this case the energy scales as 

J dal^,^-Det[vab - IpaX^X^] + ... (7.62) 

and defining dimensionless coordinates a = r^,{E)a, and dimensionless helds V* = l~^r^,{E)X'^ 
then we find 

Er^tluiE), (7.63) 

and so the Vainshtein/classicalization radius for this system is 

n{E)r.l-^E, (7.64) 

which grows with energy. The scattering amplitude is determined by the statement that 
scattering is purely elastic with the phase shift determined by the classical action 6 ~ 

S ~ Er^:{E) ~ slg with s = E'^. Since the scattering is elastic there are no ‘classicalons’ in 
this picture, nevertheless despite the perturbative non-renormalizability of the LEEFT, the 
high energy scattering amplitude is unitarized in a manner consistent with semi-classical 
arguments along the lines discussed in Sec. 2.8. 


Since the S-matrix is analytic, we may express it in terms of a dispersion relation. 
In this case, the S-matrix is polynomially bounded, since as Im{s) —>■ ±oo , —>■ 

g-|im(s)|«2/4^ In particular we can express the S-matrix in terms of a dispersion relation 
with a single subtraction 


A{s) = = 1 - - 

TT 



sin(^/^/4) 
fi{s - fl) 


(7.65) 


which defines the amplitude over the whole complex s plane, consistently with the previous 
definition. 


However as noted in [21], since the S-matrix is an entire function of 1^, there appears to 
be no logical obstruction to replacing Ig with —1^ for which the Nambu-Goto action would 
become 

S = l^g I d^a^-Det[r]ab + PgdaX%X^] . (7.66) 

This action is the analogue of the wrong sign anti-DBI model considered in Sec. 4.4.1 
for which it is known that superluminal propagation appears to arise in the LEEFT. At 
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the same time, this is the sign associated with a functioning Vainshtein mechanism. The 
S-matrix, which is now trivially 

^2i5is) ^ ^ ^ Q ^ 

g2i5(*) = ^ < 0 . (7.68) 

satisfies all of the desirable requirements of unitarity, crossing symmetry and analyticity, 
despite the apparent pathologies of the LEEFT. Once again the high energy behaviour is 
consistent with semi-classical expectations. It is precisely an analogous situation we are 
arguing is possible for Galileons with the distinction being that there we generally expect 
the scattering to be inelastic. 


We can easily see that from our language, this is an example of a non-localizable field 
theory. Indeed it is sufficient to note that in this case, the S-matrix is not polynomially 
bounded, since now as Im{s) —)> ±oo , —>■ Thus, we may still express 

the S-matrix in terms of a dispersion relation, but as in the case of all non-localizable field 
theories, we must use an entire indicator function to account for the infinite number 
of subtractions. This entire function must have the same growth at large arguments over 
the whole complex s plane as the amplitude itself. There are an infinity of such choices, 
but one consistent choice is^^ 

= sin (/i/2/4)). (7.69) 


This choice has no zeros on the real axis, and has the same exponential growth as Im(s) —)• 
Too as the S-matrix. The resulting dispersion relation is 


V (i+i (.«/4)) £ d, , (T.70) 

for a real entire function u(s) with the same order of growth. We see then that a simple 
change in the sign of /^ has a profound change of the dispersion relation due to the violation 
of polynomial boundedness. If there exists a field theory description of this system, the 
fields must necessarily be non-localizable with an order of growth a = 1. 


Following the discussion in the previous section, this model can still satisfy macro- 
locality. Reverting to non-relativistic language, this example corresponds to taking g{E) = 
E'^{1 + ^ sin (E^/^/d)). Since we are dealing with non-localizable fields of order a = 1, the 
incoming state localized at t = to should be taken to be at least as localized as a Gaussian. 

1 

*"<*) = 2ZV¥' ^ 

with L = alg and a ~ 0{1). We may view this as a smeared delta function. The response 
from this initial state is causal in terms of the modified initial state (/iin 

/ OO 

dt'G'£t - t')^£t'). (7.72) 

-OO 

^®Note that although we use the same notation, the indicator function g{E) that controls the behavior of 
the S-matrix is in general not the same as that for the Feynman propagator. 
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Here the modified initial state is 


4'm{t) = J —E^{l + -sm 

^ _(^ (2L^ -it- to)^) 
4 L5 




— 2-v/^Im 



4L^ — ill ~ “ ^o)^\ 

(4L2 - i/2)5/2) J 


(7.73) 

.(7.74) 


Since G'pj^it) = 0 for t < 0 it is clear that (l){t) will decay exponentially for t < to as 

<^out(t) ~ for t<to, (7.75) 

where 

Since the spread in energy of the incoming state is 1/7^, the region in which locality/causality 
should hold is given by the bound 

(t - to) > niE) = Ell ~ ll/L (7.77) 


and so we have in this region 

<^out(t) < exp ^-Min ^ 

and regardless of the value of L < Ig, the acausal contribution to (l)out{t) will be expo¬ 
nentially small outside the region subject to the locality bound. This is macro-causality 
without micro-causality! Despite the apparent superluminal propagation of the wrong sign 
Nambu-Goto action, Lorentz invariance and analyticity will ensure that there can be no 
macroscopic superluminal communication. 

8 Summary 

Using an explicit form for the Galileon duality transformation, we have defined a non- 
perturbative, unitary, Lorentz invariant and UV hnite quantization of a (d -|- l)’th order 
Galileon model that is dual to a free theory. We shall confirm elsewhere that the S-matrix 
is trivial [14] , establishing quantum equivalence to a free theory. By an explicit calculation 
of the Wightman functions of this theory, in two distinct duality frames, we have confirmed 
that the Kallen-Lehmann spectral densitys grow exponentially as p where r^{E) 

is the Vainshtein radius. This occurs because the resonant contribution to the spectral 
densities are dominated by V-particle states with N ~ Er^,{E). This exponential growth 
implies that this Galileon model is a Jaffe type non-localizable held theory with growth 
index a = {d + l)/(2(d -|- 2)). For such theories micro-locality and micro-causality are 
weakened to macro-locality/causality for distances larger than the locality bound 

|x|>r*(F;), (8.1) 
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where E is the typical scattering energy. This is a manifestation of the Vainshtein effect at 
the quantum level, and we conjecture that similar properties hold true for all Vainshtein 
type theories. 

Although most calculations were done using the Galileon which is dual to a free field, we 
have provided evidence that suggest all Galileon models are non-localizable field theories. 
The consequences of non-localizability for the S-matrix will be that it is no longer expected 
to be polynomially bounded, and hence Galileon scattering does not have to respect the 
Froissart bound. The consequence of this will be discussed elsewhere [14]. However, we 
have argued, following Steinmann [19], that there is no obstruction to maintaining analyt- 
icity of the retarded products, and hence analyticity of the S-matrix, despite the failure 
of polynomial boundedness. Furthermore, analyticity of the scattering amplitude alone is 
sufficient to ensure macro-causality [19]. 

Although we have explicity constructed a quantization of the dual to a free theory, the 
UV completion of an interacting Galileon remains conjectured. However, our arguments 
support the view that the UV completion cannot be a local field theory, i.e. that of a 
non-Wilsonian completion, consistent with the arguments of the ‘classicalization’ proposal 
[23-26, 45] . In particular, we note that our UV description of the special Galileon contains 
no new degrees of freedom, despite becoming strongly coupled at the scale A. In line with 
the ‘asymptotic fragility’ proposal [20, 21], our picture of the UV theory does not contain 
any UV fixed point or conformal behaviour. Rather the scale A determines the exponential 
growth properties of the Wightman functions and scattering amplitudes. We have argued 
that the superluminal group and phase velocities found in the LEEFT are generically ab¬ 
sent in the UV completion. Macro-causality will enforce luminality by modifying the high 
energy properties of the retarded propagator, which are what actually determine causal 
propagation speeds through the front velocity. 

If our conjecture is correct, then this resolves many of the technical issues with IR 
completions of Galileons such as Massive Gravity/Multi-Gravity theories [6, 12, 13]. We 
have argued that the non-locality that arises in the Galileon models is of the same nature as 
we would expect in any gravitational theory, and this fits well with the role of the Galileon as 
the helicity-zero mode of the graviton in these theories. For example, when two high energy 
particles scatter in GR at impact parameters b < Rs{E) then a black hole is expected to 
form. The production of this extended semi-classical object is tied to the mild gravitational 
non-locality of GR as proposed by Giddings and Lippert [61, 62]. When the same happens 
in Massive Gravity, as well as the Schwarzshild radius, there is the Vainshtein radius which 
describes the profile of the helicity-zero mode. The usual gravitational non-locality of GR 
is now pushed out to the larger Vainshtein/classicalization radius r^{E). In taking the 
Galileon decoupling limit [5], although the usual gravity switches off, the gravitational 
non-locality of the helicity-zero mode remains and explains why the UV completion of 
Galileons must be non-localizable. This does not undermine the classical picture of the 
Vainshtein mechanism, of relevance to cosmological phenomenology, just as black holes are 
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well-described classically. However it has profound implications for quantum fluctuations 
and the scattering of high energy Galileon quanta. 
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A A-particle phase space density 


The n-particle phase space density in 4 dimensions defined as 


can be computed as follows 


n / dfe 


, 2=1 ' 




JJ / dfci 


. 2=1 ' 




2 = 1 

g-*(*:-Er=i ki).x 
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( 4 ^ 2 )n _ (-^0 _ ’ 


(AT) 


(A.2) 

(A.3) 


Rather than computing this integral directly we can use the indirect relation for the Eu¬ 
clidean correlation functions 


1 1 


(4^2 )n y,2n 
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^oo roo 

'‘To “T (2. 

roo roo ^ 




Since we know dimensionally D„(^) = ^ then performing the /x integral 

1 1 


(4yr2)n ^2n 


1 /■“ 

= an{n-2)1^^ ds 
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B Spectral density as a Double Borel Transform 


In field theory, perturbative expansions are nearly always asymptotic and so must be 
resumed by some means, for example by Borel resummation. Our technique of defining 
Wightman functions through their spectral densities can be viewed as a double Borel 
transform. 

To see this we note that the Euclidean propagator can be written as 


roo f 1 

/ OO roo r 

ds J +^) (B.2) 

/ OO roo -I 1 

We see that the Euclidean propagator is given by a double Laplace transform of the spectral 
density. This means that the convergence properties of are always better than Ge{x). 
Specifically, if Ge{x) is given by the series 


GEix) = 


+E 


n=2 ^ ’ 


then performing the two inverse Laplace transforms we find 


^ ^ 2 (167r2)^-i(n-2)!(n-l)! 


(B.4) 


(B.5) 


This is characteristic of a double Borel transform and has a vastly improved convergence 
from the additional l/((u — 2)!(n — 1)!) at large n. In particular, if we define an operator 
via 

d(x) = (/)(x) + ^ Cn : 0(x)” :, (B.6) 

n=‘l 

then from Wick’s theorem dn = nlc^. The spectral density will converge provided only 
that Cn grows no faster than \/rd- 


C Spectral density calculation in (p, y) frame. 


A slightly different approach to the calculation of the spectral density is via a Euclidean 
calculation. The naive definition of the Euclidean two point correlation function for the 
matter field x'm. d dimensions can be written as 

d<^k 1 


{x{xi)x{x2)) = 


{2ttY 


^ik-x 


ad _ , (fe • x)^ 

X2a \y.d “ ^d+2 


(C.l) 


where x = xi —X 2 and ‘’^/^r(d/2 —1). Transforming into spherical coordinates. 


{x{xi)x{x2)) = 


(2vr) 
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where Vd -2 is the volume integral over a {d — 2) sphere and w = cos(0). Taking advantage 
of binomial expansions we may write 


(1 - dm^T = XI 


s!(n — s)! 

s=0 ^ ' 


(C.2) 


and expanding exp(ikxw) we have 




(2yr)c* _ gV \J^2aj,dJ J 

’ n=0s=0 ^ 
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dk k' 


d-\-2n—S 


E 

m=0 


{ikx) 


m rl 
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J dw - u;2)|(<^-3) . (C.3) 


Doing the w integral and then the m sum gives us 


{x{xi)x{x2)) = 


n oo n 




(27r)'^ 2^ sUn — s)\ 

^ ^ 3 n=0s=o ^ ’ 
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^2axd 


dk ^ ppQ (^s + ^-,^,s + ^-,-^k‘^x‘^^ , (C.4) 


where pFq is the regularized hypergeometric function. 

Four dimensions 

Focusing on the d = 4 case and doing the k integral we have 


. oo n / A\ 

{x{xi)x{x2)) = 7^EE 
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(27r)‘^s!(n — s)! \2tt‘^K^x^ 
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Finally, doing the s sum leaves us with 
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This series is asymptotic, as we would expect, since the real space two-point function should 
not exist even in the Euclidean. However term by term we may infer the contribution to 
the associated spectral density using the results of Appendix A to give 


. ^ ^ 47r2 “ 

P{k) = Kp) + ^ 


(27r)^ (2t^‘^A.^P n\ 

This agrees with the result used in the text. 
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